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Using  results  of  direct  numerical  stimulation  we  have  analyzed  nonlinear  energy  transfer  between  scales 
of  motion  in  isotropic  and  homogeneous  shear  turbulence  in  both  spectral  and~~physrcal  space 
representation.  In  all  cases  we  find  that  the  transfer  is  local,  occurring  between  similar  scales,  but  is 
always  moderated  by  scales  from  the  energy  containing  range  (local  transfer  through  nonlocal  triad 
interaction).  Such  local  energy  exchanges  dominate  also  subgrid-scale  energy  transfer  which  was  found 
to  be  composed  os  a  forward  and  an  inverse  transfer  components,  both  being  significant  in  dynamics  of 
large  scales.  The  spatial  structure  of  the  subgrid-scale  transfer  was  compared  with  the  structure  of  a 
number  of  physical  quantities  which  are  considered  important  in  the  dynamics  of  large  scales  giving 
usually  poor  correlations.  Moderate  correlations  were  observed  only  for  the  large  scale  energv  and  the 
Smagorinsky’s  transfer.  Finally,  a  new  theory  of  spectral  energy  dynamics  was  developed.  T  he  theor\ 
provides  a  plausible  physical,  mechanism  of  the  observed  transfer  process  and  predicts  the  energ\ 
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1  Summary  of  the  Research  Results 

The  following  goals  have  been  accomplished  in  the  course  of  this  research: 


1.  We  have  analyzed  nonlinear  energy  transfer  between  scales  of  motion  in  homogeneous 
isotropic  and  uniform  shear  turbulence  using  results  of  direct  numerical  simulations. 
The  main  conclusion  drawn  from  these  analyses  is  that  the  transfer  is  local  in  the 
spectral  space  i.e.  occurs  between  two  similar  small  scales,  but  results  from  interactions 
of  these  scales  with  much  larger  and  energetic  scales  from  the  energy  containing  range 
(local  energy  transfer  through  nonlocal  triad  interactions). 

2.  The  results  of  this  investigation  were  used  to  develop  a  new  theory  of  spectral  energy 
dynamics.  The  theory  properly  accounts  for  properties  of  interscale  energy  transfer 
observed  in  the  simulations,  predicts  the  form  of  the  energy  spectrum  in  the  dissipation 
range  in  agreement  with  experiments  and  simulations,  predicts  correctly  the  form  of 
the  energy  spectrum  in  the  inertial  range,  and  provides  a  plausible  physical  mechanism 
responsible  for  the  observed  transfer  process. 

3.  We  have  devised  a  physical  space  representation  of  the  spectral  energy  transfer  among 
scales  with  predefined  sizes  and  found  that  the  interscale  energy  transfer  is  spatially 
intermittent  and  local  in  the  physical  space. 

4.  The  physical  space  representation  of  the  spectral  energy  transfer  was  used  in  an  analysis 
of  subgrid-scale  nonlinear  interactions.  The  subgrid-scale  transfer  was  found  to  be 
composed  of  a  forward  and  an  inverse  transfer  components,  both  being  significant  in 
dynamics  of  resolved  scales.  Energy  exchanges  between  the  resolved  and  unresolved 
scales  from  the  vicinity  of  the  cutoff  wave  number  dominate  the  subgrid-scale  processes 
and  the  energetics  of  the  resolved  scales  are  unaffected  by  the  modes  with  wave  numbers 
greater  than  twice  the  cutoff  wave  number.  The  dominance  of  nonlinear  interactions 
among  the  largest  scales  in  the  subgrid-scale  energy  transfer  process  suggests  that  the 
resolved  nonlinear  term  may  serve  as  a  basis  of  a  new  approach  to  the  subgrid-scale 
modeling. 

5.  The  physical  space  representation  of  the  subgrid-scale  transfer  was  used  in  assessing 
the  dynamical  importance  of  large  scales  of  motion  (coherent  structures)  in  isotropic 
turbulence.  The  spatial  structure  of  the  exact  subgrid-scale  transfer  was  qualitatively 
compared  with  the  spatial  structure  of  a  number  of  physical  quantities  which  are  con¬ 
sidered  to  govern  the  dynamics  of  the  large  scales  of  turbulence.  It  was  found  that 
all  quantities  determined  by  the  first  derivatives  of  the  velocity  field  correlate  poorly 
with  the  transfer  which  is  largest  at  the  peripheries  of  regions  characterized  by  large 
values  of  these  quantities.  The  spatial  structure  of  the  transfer  correlates  much  better 
with  the  large  scale  energy  and  the  Smagorinsky’s  subgrid-scale  energy  transfer  which 
is  determined  by  the  second  derivatives  of  the  velocity  field.  None  of  the  considered 
quantities  was  capable  of  predicting  sign  of  the  subgrid-scale  transfer. 


2 


6.  The  interscale  energy  transfer  analysis  in  both  spectral  and  physical  representation  has 
been  performed  using  results  of  direct  numerical  simulations  of  the  Kida  flow  performed 
by  Dr.  Pelz  at  Rutgers  University.  The  results  of  this  analysis  will  serve  to  supplement 
an  ongoing  investigation  at  Rutgers  of  the  process  of  generation  of  small  scales  in  such 
a  flow. 
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5  Research  Papers 


This  section  consists  of  three  journal  papers  (one  published,  one  submitted,  and  one  in 
preparation  for  publication)  and  two  conference  proceedings,  all  in  chronological  order,  which 
provide  detailed  information  about  research  results  summarized  in  section  1. 
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Interscale  Energy  Transfer  in  Numerically 
Simulated  Homogeneous  Turbulence 

By  J.  A.  Domaradzki1,  R.  S.  Rogallo2,  and  A.  A.  Wray2 


Energy  transfer  is  investigated  for  flows  obtained  by  dirert  numerical  simulations 
of  low  Reynolds  number  homogeneous- shear  and  isotropic  turbulence  and  by  large- 
eddy  simulations  of  high  Reynolds  number  isotropic  turbulence.  The  transfer  in 
spectral  space  is  found  to  be  local  but  results  from  interaction  between  separated 
scales.  The  transfer  among  small  scales  is  highly  intermmitent  in  physical  space. 
The  measurements  suggest  an  important  correlation  between  transfer  among  small 
scales  and  the  energy  of  large  scales. 


1.  Introduction 

Using  results  of  low-Reynolds- number  direct  numerical  simulations  (DNS)  Do¬ 
maradzki  and  Rogallo  (1988,  1990)  analyzed  the  energy  transfer  in  isotropic  turbu¬ 
lence  and  concluded  that  beyond  the  energy  containing  range  the  energy  was  trans¬ 
ferred  among  scales  of  motion  similar  in  size  but  that  the  interactions  responsible 
for  this  local  energy  transfer  were  nonlocal  in  £- space.  The  same  transfer  mech¬ 
anism  was  also  found  when  the  eddy-damped  quasinormal  Markovian  (EDQNM) 
approximation  was  applied  to  high  Reynolds  number  flows  which  are  inaccessible 
to  the  DNS  technique. 

The  conclusions  concerning  the  apparent  universality  of  this  transfer  mechanism 
are  extended  in  this  work  to  homogeneous  shear  flows  and  to  high  Reynolds  number 
isotropic  flows  obtained  by  large-eddy  simulation.  We  also  devise  a  physical-space 
representation  of  the  spectral  energy  transfer  calculated  in  k  space  that  allows  us 
to  estimate  the  spatial  intermittency  of  the  energy  transfer  and  the  spatial  corre¬ 
lation  between  quantities  defined  using  only  large-scales  flow  information  and  the 
dynamically  important  energy  transfer  among  different  scales.  In  particular  this  is 
useful  in  evaluating  the  performance  of  subgrid-scale  models  formulated  in  physical 
space  e.g.  the  classical  Smagorinsky  eddy  viscosity  model. 

2.  Numerical  Velocity  Fields 

We  have  used  velocity  fields  generated  by  numerical  simulations  that  were  run 
for  sufficiently  long  times  to  fully  establish  nonlinear  interactions. 

The  velocity  field  C128U8  is  the  result  of  a  DNS  of  uniformly  sheared  homo¬ 
geneous  turbulence  performed  by  Rogers  (1986)  and  LES128  is  the  result  of  a 

1  University  of  Southern  California 
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large-eddy  simulation  of  forced  isotropic  turbulence,  at  nominally  infinite  Reynolds 
number,  performed  by  Chasnov  (1990).  Tbe  energy  spectrum  of  LES128  exhibits 
a  k~ 5/3  law  over  the  entire  range  of  simulated  wavenumbers.  The  field  K128  is 
obtained  from  a  DNS  of  isotropic  turbulence  performed  by  Rogallo  (unpublished). 
Its  use  is  motivated  primarily  by  the  fact  that  the  two  dynamically  important  pro¬ 
cesses  that  determine  the  evolution  of  the  energy  spectrum,  i.e.  viscous  dissipation 
and  nonlinear  transfer,  are  very  well  resolved.  This  resolution  is  obtained  at  the 
expense  of  lowering  the  Reynolds  number  as  compared  with  the  two  other  cases. 

3.  Basic  Quantities 

The  Navier-Stokes  equations,  in  the  Fourier  spectral  representation,  for  the  fluc¬ 
tuating  velocity  field  un  subjected  to  uniform  shear  U  =  (azj ,  0, 0)  are 

(J \  +  V&'j  »»(M)  =  (— t/2)Pnim(k)  J  it/(p,f)um(k  -  p  ,t)dp 

+2*“pr“2(k »*)  +  skl^Un^^  ~  *£u*a(k,<)  (1) 

tfcn«n  =  0  (2) 

where 

P nim(k)  =  KtiSra  -  >)  +  kt(Snm  -  knkm/k7),  (3) 

u  is  the  kinematic  viscosity,  and  the  summation  convention  is  assumed.  In  subse¬ 
quent  formulae  explicit  time  dependence  will  be  omitted. 

The  equation  for  the  energy  amplitudes  ||u(k)|2  =  |un(k)u*  (k)  is  obtained  from 

(i) 

| \ Mk)l’  =  +  T(k)  +  ~ l«(k)|J  -  ^e{Ui(k)u;(k)>  (4) 

where  the  asterisk  denotes  complex  conjugate. 

The  nonlinear  energy  transfer  is 

T (k)  =  ^  Jm  t*;(k)PnJm(k)  J  uj(p)um(k  -  p)dp  (5) 

and  the  following  two  terms  in  (4)  containing  s  describe  energy  transfer  due  to  the 
mean  shearing  deformation  of  turbulent  eddies  and  turbulent  energy  production 
by  the  mean  shear  respectively.  A  detailed  description  of  these  effects  is  given 
by  Deissler  (1961),  Fox  (1964),  and  Lumley  (1964),  and  is  summarized  in  Hinze’s 
(1975)  monograph.  Note  that  the  corresponding  equations  for  isotropic  turbulence 
are  obtained  from  (1)  by  taking  s  —  0.  In  particular,  the  nonlinear  transfer  term  (5) 
has  the  same  form  for  both  homogeneous  shear  turbulence  and  isotropic  turbulence. 


Energy  Transfer  in  Homogeneous  Turbulence 
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The  principal  quantity  of  interest  here  is  the  energy  exchange  between  a  given 
mode  k  and  all  pairs  of  modes  p  and  q  =  k  —  p  that  form  a  triangle  having  k  as 
one  of  the  legs  and  where  p  and  q  lie  in  prescribed  regions  V  and  Q  of  the  spectral 
space  respectively.  For  a  given  k,  confining  p  and  q  to  'P  and  Q  is  equivalent  to 
selecting  a  specific  set  of  triangles  from  all  of  the  possible  triangles  contributing  to 
the  energy  transfer  at  the  wavevector  k  in  (5). 

In  this  work  we  choose  V  and  Q  as  shells  in  the  wavenumber  space  k  —  |Afc  < 
|k|  <  k  +  §A k  with  a  shell  thickness  AJb.  This  choice  is  natural  for  isotropic 
turbulence  and  is  also  convenient  for  other  homogeneous  fields  as  first  suggested  by 
Batchelor  (1953). 

The  net  nonlinear  energy  transfer  to  wavenumber  band  k  is  denoted  by  T(k) 
and  the  contribution  to  this  transfer  resulting  from  nonlinear  interactions  between 
wavenumbers  in  the  band  k  and  wavenumbers  in  the  bands  p  and  q  is  denoted  by 
T(k\p,q).  According  to  this  definition 

m  =  E  £  T^> «)  =  E  w  («) 

p  «  ? 

where  the  P(k\p)  is  the  result  of  summation  of  T(k\p,q)  over  all  bands  q  and  is 
interpreted  as  the  contribution  to  the  net  energy  transfer  into  band  k  due  to  all 
interactions  involving  band  p. 

The  functions  T(Jb),  P(Jfe|p),  and  T(k\p,q)  give  progressively  more  detailed  infor¬ 
mation  about  energy  transfer  among  different  scales  of  motion  in  a  turbulent  field. 
The  method  of  computing  these  functions  is  described  by  Domaradzki  and  Rogallo 
(1990). 

4.  Analysis  of  Energy  Transfer  in  Spectral  Space 

All  of  the  contributing  terms  of  (4),  computed  for  the  field  C128U8  and  averaged 
over  spb  jrical  shells  with  thickness  A k  =  1,  are  plotted  in  figure  1.  The  calculation 
of  the  linear  transfer  skj  j|u(k)|2  suffers  from  low  accuracy  due  to  the  coarse 
resolution  of  k,  and  we  believe  that  this  term  is  close  to  zero  for  k  >  40,  contrary 
to  the  plotted  results.  Despite  this  numerical  error  a  few  important  conclusions 
can  be  drawn  from  these  results.  Nonlinear  transfer,  viscous  dissipation,  and  mean 
shear  all  make  significant  contributions  to  the  energy  balance  for  wavenumbers 
k  <  40  which  comprise  the  energy  containing  range  and  a  significant  fraction  of 
the  dissipation  range.  Energetics  of  the  smaller  eddies  (Jfc  >  40)  is  affected  only 
by  nonlinear  transfer  and  viscous  dissipation  which  are  roughly  in  balance.  Thus 
the  energetics  of  turbulence  in  about  half  of  the  spectral  domain  (k  >  40)  is  not 
affected  directly  by  the  large  scale  mean  shear. 

The  triad  structure  of  the  nonlinear  energy  transfer  term  is  illustrated  by  plotting 
P(fc|p)  in  figure  2a  as  a  function  of  k  for  p  fixed  in  a  wavenumber  band  beyond 
the  peak  of  the  energy  spectrum.  The  contributions  T(k\p,q)  to  P(ife|p),  from  all 
significant  bands  g,>*“also  included.  The  peaks  of  P(fc|p)  are  located  in  the  vicinity 
of  the  band  p  indicating  that  the  energy  transfer  is  primarily  between  comparable 
scales  of  motion.  However,  the  decomposition  into  functions  T(k\p,  q)  reveals  that 
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Figure  1.  Spectral  energy  balance  for  the  field  C128U8.  - production, 

- dissipation, - nonlinear  transfer, - linear  transfer.  The  linear  transfer 

data  has  been  smoothed. 

the  largest  contributions  to  this  local  transfer  come  from  the  interactions  involving 
a  scale  in  the  energy  containing  range  5  <  q  <  20.  Thus  for  homogeneous  shear 
flow  we  obtain  the  same  result  as  previously  reported  by  Domaradzki  and  Rogallo 
(1988, 1990)  for  isotropic  flows:  local  energy  transfer  between  two  scales  beyond  the 
energy  containing  range  results  from  nonlocal  interactions  with  scales  in  the  energy 
containing  range. 

Analysis  of  the  nonlinear  transfer  for  the  two  remaining  velocity  fields,  LES128 
and  K128,  provided  the  same  qualitative  results. 

An  attempt  was  made  to  find  a  similarity  scaling  for  the  functions  T(k\p,  q).  For 
a  given  energy  spectrum  the  following  transformation  collapses  reasonably  well  all 
curves  T(k\p,q)  for  a  band  p  beyond  the  energy  containing  range. 

T(k\p,t)  =?£(?)£(,  )T.(1^)  (7) 

The  similarity  variable  £  =  (k  —  p)/q  is  deduced  from  geometric  relations  for  a 
triad  with  legs  k,p,  and  q  and  the  scaling  factor  pE{p)E(q)  is  ad  hoc  (  but  is  found  in 
the  EDQNM  theory  for  power-law  spectra  in  the  disparate-scale  limit ).  In  figure  2b 
we  show  the  result  of  scaling  (7)  applied  to  the  measured  functions  T(k\p,  q)  of 
figure  2a.  Interestingly,  the  transfer  scales  with  the  energy  E(q)  of  the  large  eddies 
rather  than  with  their  rate-of-strain  qE(q)lt*  which  is  the  scaling  postulated  by  a 
number  of  classical  closure  hypotheses  (Monin  and  Yaglom,  1975).  We  have  not 
been  able  to  propose  a  convincing  dynamical  model  of  transfer  processes  which 
would  provide  scaling  (7). 
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Figure  2.  Detailed  triad  contributions  to  energy  transfer  for  case  C128U8:  (a) 
unsealed,  (b)  scaled  by  (7).  The  transfer  spectra  T(k\p,q)  are  shown  for  band 
40  <  p  <  45,  and  all  bands  q  that  make  a  significant  contribution  to  P(k\p). 

H - 0  <  q  <  5,  a - 5  <  q  <  10,  □ .  10  <  q  <  15,  a .  15  <  q  <  20, 

H - 20  <  q  <  25, - P(Jfc|p). 
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5.  Physical  Space  Representation  of  Spectral  Energy  Transfer 

Let  us  denote  by  N%Q(k)  the  contribution  to  the  integral  (the  nonlinear  term) 
in  (1)  from  only  those  interactions  between  modes  p  and  q  =  k  —  p  such  that 
each  of  them  is  confined  to  one  of  the  two  prescribed  wavenumber  bands  V  and  Q. 
This  quantity  is  computed  using  the  method  described  by  Domaradzki  and  Rogallo 
(1990).  Its  Fourier  transform  to  physical  space,  N^Q(x)  say,  gives  the  contribution 
to  the  rate  of  change  of  velocity  in  physical  space  un(x,<)  caused  by  the  nonlinear 
interactions  involving  two  scales  from  the  respective  wavenumber  bands  V  and  Q 
in  the  spectral  space.  Note  that  these  interactions  influence  all  modes  k  that  can 
form  a  triangle  with  modes  such  that  one  is  in  "P  and  the  other  in  Q.  Consider  next 
a  velocity  field  truncated  to  a  prescribed  wavenumber  band  K,  i.e. 


..X/kx  _  f  wn(k),  if  k eK 
n  .  \  0,  otherwise. 


(8) 


The  Fourier  transform  of  (8)  to  physical  space,  uj(x)  say,  represents  the  contri¬ 
bution  in  physical  space  that  scales  from  band  K  make  to  the  total  velocity.  The 
contracted  product  of  these  two  physical  space  quantities 


TKTe(x)  =  ui(x)NZa(x)  (9) 

gives  a  physical  space  representation  of  the  energy  transfer  to/from  modes  in  the 
fc-band  due  to  their  nonlinear  interactions  with  modes  in  the  p-  and  g-bands. 

An  interesting  case  is  obtained  by  dividing  wavenumber  space  into  two  disjoint 
regions  K  (k  <  kc)  and  V  (k  >  kc).  The  quantity 

r5G5(x|fcc)  =  tkvv{x)  +  tkkv(x)  (10) 

provides  a  physical  space  representation  of  the  rate  of  change  of  energy  of  large 
scales  k  <  kc  due  to  nonlinear  interactions  involving  small  scales  k  >  kc.  This  is 
precisely  the  energy  transfer  process  which  is  the  suo  ject  of  subgrid-scale  modelling. 

We  have  computed  transfer  functions  (9)  and  (10)  for  various  wavenumber  bands 
of  the  field  K128.  The  low  wavenumber  band  Q  is  chosen  to  cover  the  entire 
energy  containing  range  0  <  q  <  10.  Figure  3a  shows  one  plane  from  the  full 
transfer  (9)  representing  in  physical  space  the  energy  transfer  to  eddies  in  the  band 
23  <  A:  <  28  caused  by  their  interactions  with  eddies  in  the  bands  20  <  p  <  25  and 
0  <  q  <  10.  The  transfer  function  is  spatially  intermittent  and  is  predominantly 
positive  indicating  a  flow  of  energy  from  the  larger  scales  p  to  the  smaller  scales 
k.  In  figure  4b  we  plot  the  same  function  for  17  <  k  <  22.  The  transfer  is  now 
predominantly  negative  as  expected  and  occurs  at  roughly  the  same  locations  as 
the  transfer  of  figure  3a.  We  thus  conclude  that  the  local  energy  transfer  between 
similar  wavenumber  modes  in  spectral  space  is  intermittent  in  physical  space. 

We  have  attempted  to  correlate  this  spatial  distribution  of  energy  transfer  with 
a  number  of  simpler  quantities  (rate-of-strain,  dissipation,  energy,  etc.)  calculated 
from  the  velocity  field  truncated  to  contain  only  either  large  or  small  scales.  In 
figures  4a  and  4b  we  show  the  physical-space  distribution  of  energy  for  the  velocity 
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field  truncated  to  23  <  k  <  28  and  0  <  k  <  10,  respectively.  Both  energy  fields 
correlate  very  well  with  the  energy  transfer  among  small  scales  shown  in  figure  3. 
Correlation  of  other  calculated  quantities  with  the  energy  transfer,  notably  the 
square  of  the  rate-of-strain  tensor,  was  generally  much  worse.  Therefore  we  conclude 
that  the  energy  transfer  among  small  scales  occurs  mostly  at  those  physical  locations 
that  contain  large  amounts  of  turbulent  energy  rather  than  at  the  locations  of 
high  strain  rate.  This  correlation  is  the  physical  space  counterpart  of  the  observed 
importance  of  the  nonlocal  triads  in  the  energy  transfer  process  in  spectral  space. 

We  have  used  formula  (10)  to  calculate  subgrid-scale  (SGS)  energy  transfer  for 
the  field  K128  with  the  cutoff  wavenumber  ke  —  10.  The  full  SGS  transfer  field, 
plotted  in  figure  5a  for  a  typical  plane,  is  characterized  by  the  presence  of  both 
negative  and  positive  regions.  These  indicate  energy  flux  from  and  to  the  large 
scales  respectively  due  to  subgrid-scale  interactions.  The  classical  Smagorinsky 
model  (Smagorinsky,  1963)  for  this  transfer,  based  on  the  velocity  field  truncated 
to  the  large  scales  0  <  k  <  10,  is  plotted  in  figure  5b.  Note  that  the  model  captures 
properly  the  locations  of  the  regions  where  the  transfer  is  most  intense  but  fails 
completely  to  predict  the  inverse  energy  transfer  from  small  to  large  scales. 

6.  Conclusions 

Using  results  of  direct  numerical  simulations  of  homogeneous  shear  turbulence  we 
have  shown  that  the  nonlinear  energy  transfer  in  spectral  space  beyond  the  energy 
containing  range  has  the  same  character  as  reported  previously  for  isotropic  turbu¬ 
lence:  local  energy  transfer  caused  by  nonlocal  triad  interactions.  The  same  conclu¬ 
sion  was  reached  for  velocity  fields  obtained  in  large-eddy  simulations  of  isotropic 
turbulence  at  high  Reynolds  numbers. 

An  ad  hoc  scaling  roughly  collapses  the  transfer  T(k\p,q)  to  a  self-similar  form. 
This  scaling  implies  an  important  role  which  the  energetic  scales  play  in  the  energy 
transfer  among  small  scales,  but  the  process  does  not  appear  to  be  simply  straining 
of  the  small  scales  by  the  large  ones. 

We  have  devised  a  physical  space  representation  of  the  energy  transfer  processes 
among  scales  of  motion  belonging  to  three  distinct  wavenumber  bands  in  spectral 
space  and  conclude  from  it  that  the  energy  transfer  among  small  scales  is  highly 
intermittent  in  physical  space.  Furthermore,  regions  of  significant  transfer  appear 
to  correlate  better  with  regions  of  significant  large-scale  energy  than  with  those  of 
significant  large-scale  strain  rate. 

As  a  particular  case  we  have  calculated  the  subgrid-scale  energy  transfer  in 
isotropic  turbulence.  This  SGS  transfer  exhibits  regions  of  energy  drain  from  large 
to  small  scales  as  well  as  significant  regions  of  reversed  energy  transfer  from  small 
to  large  scales.  The  Smagorinsky  eddy  viscosity  model  captures  the  locations  of 
the  most  intense  transfer  but  predicts  that  it  is  always  from  large  to  small  scales, 
contrary  to  the  measurements  from  direct  calculations. 
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Figure  5.  Subgrid  scale  energy  transfer  in  physical  space  Tsgs  (x|Jbe)  for  Jfec  =  10: 
(a)  measured,  (b)  computed  using  the  Smagorinsky  eddy  viscosity  model  for  the 
velocity  field  truncated  to  wavenumber  band  0  <  Jfe  <  10. 
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Abstract-Direct  numerical  simulations  of  homogeneous  turbu¬ 
lence  are  used  to  analyze  energy  transfer  among  scales  of  motion  in 
spectTal  space.  A  physical  space  representation  of  such  a  spectral 
energy  transfer  is  devised  and  applied  to  the  analysis  of  an  eddy 
viscosity  with  a  sharp  spectral  cut-off. 

1  Introduction 


where  NnPa  is  (3)  calculated  with  one  of  the  contributing  velocity 
fields  truncated  to  V  and  the  other  to  Q.  Details  of  such  calculations 
are  provided  in  [2].  For  homogeneous  turbulence  the  regions  V  and 
Q  are  usually  chosen  as  spherical  wavenumber  bands.  Similarly  trun¬ 
cating  velocity  u*  (k)  in  (4)  to  a  spherical  shell  X  results  in  a  quantity 
TKva(k)  which,  after  averaging  over  X,  is  interpreted  as  the  energy 
transfer  to  the  band  X  resulting  from  nonlinear  interactions  of  scales 
in  X  with  scales  in  V  and  £■ 


Statistically  homogeneous  turbulent  flows  are  conveniently  repre¬ 
sented  in  spectral  (Fourier)  space.  In  such  a  representation  dynami¬ 
cally  important,  elementary  nonlinear  interactions  involve  three  dis¬ 
tinct  modes  with  their  wavenumbers  forming  a  closed  triad.  Under¬ 
standing  these  interactions  is  of  paramount  importance  in  the  theory 
of  turbulence  since  essentially  all  turbulence  closures  rely  on  assump¬ 
tions  about  the  nature  of  the  nonlinear  interactions.  Recently,  using 
results  of  direct  numerical  simulations  (DNS)  Domaradzlti  and  Ro- 
gallo  [1]  [2J  analyzed  the  energy  transfer  in  homogeneous  turbulence. 
They  concluded  that  beyond  the  energy  containing  range  the  energy 
was  transferred  among  scales  of  motion  similar  in  size  but  that  the 
interactions  responsible  for  this  local  energy  transfer  were  nonlocal 
in  L space.  The  importance  of  such  nonlocal  triadic  interactions  in 
the  evolution  of  turbulent  flows  has  been  confirmed  by  Yeung  and 
Brasseur  [3j  who  also  provided  analytical  arguments  [4)  supporting 
conclusions  drawn  from  DNS. 

Despite  the  usefullness  of  spectral  representation  as  a  theoretical 
and  numerical  tool  in  turbulence  research,  various  quantities  (veloc¬ 
ity,  energy,  varticity,  etc.)  in  the  physical  space  often  provide  a  more 
natural  description  of  turbulent  flows.  Thus  it  is  of  interest  to  have 
the  physical  space  representation  of  the  nonlinear  transfer  processes 
that  dominate  the  spectral  space  dynamics.  One  such  representation 
has  been  proposed  by  Domaradzlti  et  aL  [5].  In  this  paper  we  discuss 
other  possible  ways  of  representing  detailed  spectral  energy  transfer 
in  the  physical  space. 

2  Interscale  Energy  Transfer  in  Spectral  Space 

The  equation  for  the  energy  amplitudes  l|u(k)|J  =  Iun(k)u;(k) 

is: 


3  Interscale  Energy  Transfer  in  Physical  Space 


Inverse  Fourier  transform,  signified  by  tilde,  of  N„(k)  is  the  sum 
of  the  convection  and  pressure  terms  in  the  Navier-Stokes  equation 
in  the  physical  space  coordinates 


N„(x)  =  -tii(x) 


dun(x) 

dzi 


opjx) 

dz„ 


(5) 


Similarly,  using  we  can  define  its  physical  space  coun¬ 

terpart  JV^a(x)  as  well  as  Nr^”e(x)  which  is  the  inverse  Fourier 
transform  of  Nn’,c(k)  truncated  to  the  band  X.  N%a{x)  can  be 
interpreted  as  the  contribution  to  the  rate  of  change  of  the  velocity 
field  {^(x)  at  a  point  x  made  by  the  nonlinear  interactions  involving 
modes  from  the  bands  V  and  Q.  Note  that  these  interactions  influ¬ 
ence  all  modes  k  which  c»n  form  a  triangle  with  two  other  inodes 
such  that  one  is  in  V  and  the  other  in  Q.  jVjf^x)  represents  a  con¬ 
tribution  to  the  rate  of  change  of  u„(x)  which  is  made  by  all  inodes 
from  X  interacting  nonlinearly  with  modes  in  P  and  Q- 

The  rate  of  change  of  the  turbulent  energy  e(x)  =  ^u»(x)tb«(x) 
at  a  point  x  caused  by  the  nonlinear  interactions  is 


£&«*.(*)&(*)•  (6) 

Our  goal  is  to  decompose  (6)  into  contributions  from  the  interac¬ 
tions  among  modes  from  predefined  wavenumber  bands  *  and  Q 
i.e.  to  find  a  physical  space  counterpart  of  k)  which  itself  is 
the  result  of  such  a  decomposition  of  the  transfer  T(k)  p  erf  armed  in 
the  spectral  space.  Despite  uniqueness  of  such  a  decomposition  in  the 
spectral  representation*  the  procedure  is  ambiguous  in  the  physical 
space.  Possible  definitions  are: 


|i|u(k)|J  =  -2Vtli|u(k)|I  +  r(k)  (i) 

where  u„(k)  i«  the  velocity  field  in  spectral  space,  with  the  ex¬ 
plicit  dependence  on  time  omitted,  the  asterisk  denotes  complex  con¬ 
jugate,  v  is  the  kinematic  viscosity,  and  T(k)  is  the  nonlinear  energy 
transfer 

T(k)  =  Jl*«(k)JV4k)).  (2) 

In  the  last  equation  Nn(k)  is  the  nonlinear  term  in  the  Navier- 
Stokes  equation 

N*(k)  =  (-i/2)JWk)  j  <*WpK.(k  -  p) ,  (3) 

where  tensor  P«im(k)  account!  for  the  pressure  and  incompress¬ 

ibility  effects.  The  summation  convention  is  assumed  throughout. 

Detailed  energy  transfer  to/from  mode  k  caused  by  its  interac¬ 
tions  with  wavenumbers  p  in  a  prescribed  region  “P  of  the  wavenum¬ 
ber  space  and  q  =  k  —  p  in  another  region  Q  is 

T*a(k)  =  fleK(k)N/a(k))  (4) 


TiCTC(x)  =  ^(x)Nn”a(x) ,  (7) 

fsKTQ(x)  =  «*(x)l9**8(x) ,  (8) 

f~a(x)  =  u»(x)Nf’a(x) ,  (9) 

where  uj(x)  is  the  inverse  Fourier  transform  of  Un(k)  truncated 

to  the  band  X.  —jcrQiiA 

Function  is  a  straightforward  counterpart  of  “  (*»’ 

with  a  product  of  fi*  and  N*a  taken  in  the  physical  rather  than  ® 
the  spectral  space.  However,  since  u.K(k)  vanishes  outride  X.tM 
multiplication  in  the  spectral  space  implicitly  truncates  Nn  (  ) 
the  band  so  that  TCPa(k)  expresses  transfer  to  the  n»desm 
X  only.  In  f,CTa(x)  the  effect  of  nonlinear  transfer  to  modes  onto 
X  is  present  in  the  tenn  JVjf  °.  . 

An  explicit  truncation  of  N.,a(k)  to  X  and  multiplication  oj 
uf  seems  to  rectify  this  problem  resulting  in  (8).  The  drawback 
this  definition  L-  that  it  does  not  satisfy  a  natural  condition: 

£  r"a(x)  =  2£i.  (i») 
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which  is  satisfied  by  both  (7)  and  (^9). 

Function  ?"e(x)  may  be  interpreted  as  a  fraction  of  the  rate 
of  change  of  the  total  energy  e(x)  due  to  variation  of  modes  in  1C  as 
they  are  affected  by  nonlinear  interactions  with  modes  from  V  and 
2- 

Thus  none  of  the  above  definitions  is  an  exact  counterpart  of  the 
spectral  transfer  T**®(k)  but  (9)  is  the  most  appealing  candidate. 

An  interesting  special  case  is  obtained  by  dividing  a  wavenumber 
space  into  two  disjoint  regions  K.  (k  <  kc)  and  V  {k  >  kc).  Quantity 

Tscstxii,)  =  +  TKK*(x)  (11) 

provides  a  physical  space  representation  of  the  rate  of  change  of 
energy  of  large  scales  (k  <  kc)  due  to  their  nonlinear  interactions 
through  wavenumber  triads  which  have  at  least  one  of  the  legs  in  the 
region  P.  This  is  precisely  the  energy  transfer  process  which  is  the 
subject  of  the  subgrid-scale  modeling. 

We  have  computed  transfer  functions  (7)  and  (11)  for  the  statisti¬ 
cally  isotropic  velocity  field  obtained  in  direct  numerical  simulations 
performed  with  a  resolution  of  12S3  modes  (maximum  wavenumber 
k  =  64).  The  low  wavenumber  band  Q  remains  always  fixed  and  is 
chosen  to  cover  the  entire  energy  containing  range  (0  <  q  <  10).  Fig¬ 
ure  1  shows  one  plane  from  the  full  transfer  (7)  representing  in  the 
physical  space  the  energy  transfer  to  eddies  in  the  band  23  <  i  <  28 
caused  by  their  interactions  with  eddies  in  the  bands  20  <  p  <  25 
and  0  <  q  <  10.  The  transfer  function  is  spatially  intermittent  and 
is  predominantly  positive,  indicating  a  flow  of  energy  from  the  larger 
scales  p  to  the  smaller  scales  k. 

We  have  attempted  to  correlate  this  physical  energy  transfer  with 
a  number  of  simpler  quantities  (rate-of-strain,  dissipation,  energy, 
etc.)  calculated  from  the  velocity  field  truncated  in  such  a  way  as  to 
contain  only  either  large  or  small  scales.  We  found  that  the  energy  of 
the  velocity  field  truncated  to  large  scales  0  <  i  <  10  correlates  very 
well  with  the  energy  transfer  among  small  scales  shown  in  figure  1. 
Correlation  of  other  calculated  quantities  with  the  energy  transfer, 
notably  the  square  of  the  rate-af-strain  tensor,  was  generally  much 
worse.  Therefore  we  conclude  that  the  energy  transfer  among  small 
scales  occurs  mostly  at  those  physical  locations  which  contain  large 
amounts  of  turbulent  energy  rather  than  at  the  locations  of  high 
strain  rate,  an  unexpected  result.  Indeed,  until  this  paradox  is  re¬ 
solved,  we  can  not  be  confident  that  the  particular  measure  of  energy 
transfer  that  we  have  used  is  the  appropriate  one. 

We  have  used  formula  (11)  to  calculate  subgrid-scale  (SGS)  en¬ 
ergy  transfer  for  the  same  field  with  the  cutoff  wavenumber  be  =  10. 
A  plane  from  the  full  SGS  transfer  field  is  plotted  in  figure  2.  The 
transfer  is  characterized  by  the  presence  of  both  negative  and  positive 
regions.  These  indicate  energy  flux  from  and  to  the  Large  scales  re¬ 
spectively  due  to  their  interactions  with  the  smaller  scales.  Standard 
subgrid-scale  eddy  viscosity  models  predict  transfer  in  one  direction 
only,  from  large  to  small  scales. 

4  Conclusions 

We  have  devised  a  physical  space  representation  of  the  energy 
transfer  processes  among  scales  of  motion  belonging  to  three  distinct 
wavenumber  bends  in  the  spectral  space  and  conclude  from  it  that 
the  energy  transfer  among  small  scales  is  highly  intermittent  in  the 
physical  space  and  correlates  well  with  regions  of  significant  large- 
scale  energy. 

As  a  particular  case  we  have  calculated  a  subgrid-scale  energy 
transfer  in  isotropic  turbulence.  The  SGS  transfer  exhibits  regions 
of  energy  drain  from  large  to  small  scales  ss  well  as  significant  regions 
of  reversed  energy  transfer  from  small  to  large  scales.  Classical  eddy 
viscosity  models  assumt  that  transfer  is  always  from  large  to 
scales,  contrary  to  the  results  of  direct  calculations. 
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Fig.l.  Energy  transfer  in  physical  space  TKVa  for  23  <  k  < 
28,20  <p<  25,0  <?<  10. 


Fig.2.  Subgrid  scale  energy  transfer  in  physical  space  Tsoi(x|t*) 
for  kc  =  10. 
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Detailed  transfer  functions  T(k\p,q),  which  express  turbulent  energy  transfer  rate  to  modes 
k  caused  by  their  nonlinear  interactions  with  modes  p  and  q,  are  analyzed  using  results 
of  direct  numerical  simulations  of  homogeneous  turbulence.  A  previously  found 
phenomenological  scaling  for  the  functions  T{k\p,q),  which  brings  them  into  a  self-similar 
form,  is  used  to  deduce  the  form  of  the  energy  spectrum  in  the  dissipation  range 
proportional  to  k~z  exp(  —ak)  and  the  transfer  spectrum  proportional  to  exp(  —ak).  A 
physical  mechanism  of  the  energy  transfer  process  consistent  with  the  self-similarity  scaling  is 
proposed. 


I.  INTRODUCTION 

Navier-Stokes  equations  for  an  incompressible,  homo¬ 
geneous  turbulent  flow  have  the  following  form  in  spectral 
(Fourier)  representation  (see,  e.g.,  Lesieur,1  pp.  92-94): 

d 

^«,(k)  =  -vfc2u,(k)+iV„(k).  (1) 

Here,  u„(k)  is  the  velocity  field  in  spectral  space,  with  the 
explicit  dependence  on  time  omitted,  v  is  the  kinematic 
viscosity,  and  .V„(k)  is  the  nonlinear  term 

JV„(k)  =  (~jpn,„(k)  J  dlp  H,(p)«m(k-p),  (2) 

where  tensor  Pnlm( k)  accounts  for  the  pressure  and  incom¬ 
pressibility  effects.  The  summation  convention  is  assumed 
throughout. 

The  equation  for  the  energy  amplitudes  j|u(k)|2 
=  ju,(k)uJ(k),  where  the  asterisk  denotes  a  complex  con¬ 
jugate,  is  obtained  from  ( 1 )  and  has  the  following  form: 

<9  1  1 

jtj \u(k)\2=-lvk>- |a(k)|24-r(k).  (3) 

In  the  last  equation  Hit)  is  the  nonlinear  energy  transfer 

r(lc)  =Re[M*(k)AJ,(k)  ],  (4) 

expressing  the  rate  of  change  of  energy  of  the  mode  k 
caused  by  its  nonlinear  interactions  with  all  other  modes  in 
the  system. 

Detailed  energy  transfer  to/from  mode  k  caused  by  its 
interactions  with  wave  numbers  p  in  a  prescribed  region  SP 
of  the  wave-number  space  and  q=k— p  in  another  region 
SI  is 

r^(k)=ReK(k)JVf2(k)],  (5) 

where  Wj^is  given  by  formula  (2),  calculated  with  one  of 
the  contributing  velocity  fields  truncated  to  SP  and  the 
other  to  S2 .  Details  of  such  calculations  are  provided  by 
Domaradzki  and  Rogallo.2*3  For  homogeneous  turbulence 
the  regions  S?  and  SI  are  usually  chosen  as  spherical  wave¬ 


number  bands  centered  at  wave  numbers  p  and  q,  respec¬ 
tively.  Similarly,  truncating  velocity  u*{ k)  in  (5)  to  a 
spherical  shell  JF,  centered  at  a  wave  number  k,  results  in 
a  quantity  7^r^i>(k),  which  is  nonzero  only  for  those 
modes  k  that  belong  to  the  region  JC .  Summing  this  last 
quantity  over  all  modes  in  JF  results  in  a  function  denoted 
as  T{k\p,q),  which  is  interpreted  as  the  energy  transfer  to 
the  band  JS  resulting  from  nonlinear  interactions  of  scales 
in  JF  with  scales  in  SP  and  S>. 

Total  nonlinear  energy  transfer  T{k)  to  a  wave- 
number  band  k  (i.e.,  the  region  JT)  is  obtained  by  sum¬ 
ming  contributions  T{k\p,q)  from  all  possible  bands p  and 
V 

rot)—  S  I  nk\p,q)=^p{k\P).  (6) 

pi  p 

Here,  the  function  P{k\p)  is  a  result  of  summation  of 
T(k\p,q)  over  all  q  bands  and  is  interpreted  as  the  energy 
transfer  between  wave-number  bands  k  and  p.  Also  note 
that  T(k)  can  be  obtained  from  (4)  by  summing  T{ k) 
over  all  modes  in  the  wave-number  shell  k. 

Functions  T{k),  P(k\p),  and  T(k\p,q)  give  progres¬ 
sively  more  detailed  information  about  energy  transfer 
among  different  scales  of  motion  in  the  turbulent  field. 

For  low  Reynolds  number  homogeneous,  isotropic, 
and  shear  turbulence,  the  same  qualitative  result  for  the 
function  T{k\p,q)  was  reported  by  Domaradzki  and 
Rogallo2'3  and  Yeung  and  Brasseur*  (isotropic  turbu¬ 
lence),  and  by  Domaradzki  et  al 5  (shear  turbulence):  lo¬ 
cal  energy  transfer  between  two  scales  k  and  p  outside  the 
energy  containing  range  caused  by  nonlocal  interactions 
with  the  third  scale  q  in  the  energy  containing  range.  For 
high  Reynolds  number  flows,  similar  behavior  of  the  func¬ 
tion  T(k\p,q)  was  observed  by  Ohkitani  and  Kida,6  using 
results  of  direct  numerical  simulations  of  a  high  symmetry 
flow,  and  by  Domaradzki  et  aL ,5  using  velocity  fields  ob¬ 
tained  in  large  eddy  simulations  of  Chasnov.7  Moreover, 
the  function  T(k\p,q)  computed  by  Domaradzki  and 
Rogallo3  and  Ohkitani  and  Kida6  in  the  framework  of  the 
analytical  theories  of  turbulence  exhibits  the  same  features 
as  observed  in  the  simulations. 
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Even  though  there  is  little  disagreement  concerning 
properties  of  the  function  T(k\p.q).  there  is  a  serious  dis¬ 
agreement  concerning  physical  interpretation  and  signifi¬ 
cance  of  the  observed  form  of  T{k\p,q).  in  particular,  for 
the  energy  transfer  process  in  the  inertial  range  of  turbu¬ 
lence.  Using  an  asymptotic  analysis  of  mad  interactions, 
Brasseur8  argued  that  the  dynamical  couplings  between 
large  and  small  scales  strengthen  with  increasing  Reynolds 
number,  and  this  effect  casts  doubt  on  the  validity  of  the 
classical  assumptions  of  local  isotropy  and,  consequently, 
may  require  modifications  in  Kolmogorov’s  arguments, 
leading  to  the  universal  form  of  the  inertial  range  spec¬ 
trum.  On  the  other  hand,  Waleffe9  dismisses  physical  sig¬ 
nificance  of  such  interactions  entirely,  asserting  that  "... 
the  nonlocal  interactions  with  local  transfer  character  of 
triadic  interactions  is  not  property  of  turbulence  physics, 
but  rather  a  general  feature  of  the  Fourier  representation.” 
Others  (Zhou  and  Rogallo10)  are  inclined  to  treat  as  phys¬ 
ically  interpretable  quantities  only  certain  integrals  of 
T(k\p,q),  e.g.,  net  energy  flux  across  the  spectrum,  and 
not  the  function  T{k\p,q)  itself. 

To  some  extent,  these  controversies  are  caused  by  a 
fairly  qualitative  character  of  the  above  referenced  analyses 
and  our  poor  understanding  of  the  relation  between  dual 
pictures  of  turbulence,  one  using  a  physical  space  and  the 
other  a  spectral  (Fourier)  space  representation.  In  this 
paper  we  draw  several  quantitative  conclusions  about  the 
energy  and  the  transfer  spectra  from  the  observed  form  of 
the  function  T(k\p,q)  and  propose  a  particular  mecha¬ 
nism  of  interactions  between  scales  in  the  physical  space, 
which  is  consistent  with  the  observed  behavior  of  the  func¬ 
tion  T(k\p,q)  in  the  spectral  space. 

II.  SELF-SIMILARITY  SCALING  FOR  THE  FUNCTION 

T(k\p,q) 

In  Ref.  5  three  different  velocity  fields  generated  by 
numerical  simulations  were  considered.  The  velocity  field 
C128U8  was  the  result  of  a  direct  numerical  simulation  of 
uniformly  sheared  homogeneous  turbulence  performed  by 
Rogers  et  al.u  and  LES128  was  the  result  of  a  large-eddy 
simulation  of  forced  isotropic  turbulence,  at  a  nominally 
infinite  Reynolds  number,  performed  by  Chasnov.7  The 
energy  spectrum  of  LES128  exhibits  the  Kolmogorov 
k~i/i  law  over  the  entire  range  of  simulated  wave  num¬ 
bers.  The  field  K128  was  obtained  from  a  direct  numerical 
simulation  of  isotropic  turbulence  performed  by  Rogallo. 12 

For  all  three  fields  the  functions  T{k\p,q)  were  com¬ 
puted.  In  Fig.  1(a)  we  show  the  functions  T(k\p,q)  com¬ 
puted  for  the  field  C128U8  for  p  fixed  outside  the  energy- 
containing  range.  The  triad  structure  of  the  nonlinear 
energy  transfer  term  is  best  understood  by  considering  in 
this  figure  P(k\p)  as  a  function  of  k  (for  p  fixed)  and 
decomposing  P{k\p)  into  functions  T(k\p,q)  for  all  pos¬ 
sible  wave-number  bands  q.  Peaks  of  P{k\p)  are  located  in 
the  vicinity  of  the  prescribed  p  band,  indicating  that  the 
energy  transfer  is  most  effective  among  comparable  scales 
of  motion.  However,  the  decomposition  into  functions 
T(k\p,q)  reveals  that  the  largest  contributions  to  this  local 
transfer  come  from  the  triad  interactions  with  the  third 


(a) 


4CC 
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FIG.  1.  Detailed  energy  transfer  functions  T{k\p,q)  for  field  C128U8; 
(a)  unsealed,  (b)  scaled  using  Eq.  (7).  The  transfer  spectra  T{k\p,q)  are 
shown  for  band  40</><45,  and  all  bands  q  that  make  a  significant  con¬ 
tribution  to  P{k\p)  [the  solid  line  in  (a)).  The  vertical  dashed  lines 
delineate  the  band  p  and  different  bands  q  for  which  transfer  curves 
T(k\p,q)  are  plotted:  dashed  line  and  plus  symbols.  0  < 9  < 5;  double- 
dashed  line  and  triangles,  5  <  q  <  10;  dotted  line  and  squares.  10  <  q  <  1 5; 
dotted  line  and  triangles,  15<9<20;  double-dashed  line  and  plus  sym¬ 
bols,  20<9<25.  In  (a)  peaks  of  the  transfer  curves  broaden  for  increas¬ 
ing  q. 


interacting  scale  in  the  energy-containing  range  (here,  5 
<q  <  20).  Thus  we  observe  local  energy-  transfer  between 
two  scales  outside  the  energy-containing  range  caused  by 
nonlocal  interactions  with  the  third  scale  in  the  energy- 
containing  range. 

In  Ref.  5  it  was  found  that  for  a  given  energy  spectrum 
the  following  transformation  collapsed  reasonably  well  all 
curves  T(k\p,q)  for  all  bands  p  outside  the  energy- 
containing  range: 


‘\  q  j  ~ pE{p)E(q) 


(7) 


The  similarity  variable  £=(fc—  p)/q  is  deduced  from 
geometric  relations  for  triads  with  legs  k  p,  and  q.  For  a 
fixed  p  and  q,  with  p>q,  the  range  of  wave  numbers  k  that 
can  form  a  triangle  with  those  wave  numbers  is  p—q<k 
<P+q •  Therefore,  for  any  such  pair  of  fixed  wave  numbers 
p  and  q  the  variable  £  will  always  be  in  the  range  —  1  <  £ 
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<  1'.  The  scaling  factor  pE(piE{q'  is  phenomenological. 
The  result  of  scaling  i  "  i  applied  to  the  functions  T'  k  p.q) 
.is  plotted  in  Fig.  l(bi.  The  quality  of  the  scaling  was  com¬ 
parable  for  all  three  fields  considered  and  it  was  the  only 
scaling,  among  those  tned.  that  had  any  success  in  collaps¬ 
ing  the  numerical  data.  We  must  emphasize  that  the  above 
relation  applies  to  a  particular  combination  of  two  fixed 
wave  numbers,  q  being  in  the  energy-containing  range  and 
p  >  q  outside  this  range.  Therefore,  even  though  the  func¬ 
tion  T(k\p,q)  is  symmetric  in  p  and  q  the  scaling  is  not, 
since  this  particular  ordering  of  the  wave  numbers  p  and  q 
is  assumed  in  (7). 

No  explanation  for  the  scaling  (7)  exists  and  at  this 
stage  it  must  be  treated  as  a  purely  phenomenological  re¬ 
sult  akin  to  experimental  fitting  procedures.  Scatter  in  the 
data  in  Fig.  1(b)  is  caused  partially  by  this  approximate 
nature  of  the  scaling  and  partially  by  the  fact  that  in  our 
numerical  procedure  we  are  dealing  with  wave-number 
bands  of  finite  thickness  (A  k  —  l,Ap=  Aq  =  5)  and  not  with 
sharply  defined  wave  numbers  k.  p,  and  q.  Despite  these 
shortcomings  of  the  scaling  relation  (7)  in  what  follows  we 
will  assume  that  it  holds  for  the  functions  T{k  \p,q)  and  we 
will  explore  its  consequences  and  a  possible  physical  inter¬ 
pretation. 

It  is  important  to  note  that  according  to  (7)  the  trans¬ 
fer  scales  formally  with  the  energy  E(q)  of  the  large  eddies 
rather  than  with  their  rate  of  strain  qE(q)l/z,  which  is  the 
scaling  postulated  by  a  number  of  classical  closure 
hypotheses. 13  At  the  present  time  the  dependence  of  the 
energy  transfer  process  among  small  scales  on  the  energy  of 
large  scales  observed  in  the  simulations  lacks  a  clear  phys¬ 
ical  explanation.  Direct  coupling  between  large  and  small 
scales  is  postulated  in  the  linear  theories  of  the  viscous- 
convection  range  of  Batchelor14  and  Kraichnan15  and  the 
far  dissipation  range  of  Novikov16  and  Saffrnann,17  but 
these  theories  assume  that  the  small  scales  are  affected  by 
the  strain  of  the  large  scales,  not  their  energy.  For  turbu¬ 
lence  in  the  inertial  range  the  classical  theories  predict  that 
the  energy  transfer  is  dominated  by  the  local  interactions. 
For  instance,  theories  of  Obukhov  and  of  Heisenberg  ( see 
Monin  and  Yaglom13)  express  a  spectral  energy  flux 
through  a  wave  number  k  as  a  product  of  two  wave- 
number  integrals  involving  the  energy  spectrum  £(k)  with 
the  largest  contributions  to  both  integrals  coming  from  the 
wave  numbers  in  the  vicinity  of  k.  Thus  the  nonlocal  in¬ 
teractions  discussed  here  are  considered  to  be  of  little  im¬ 
portance  in  those  theories.  Similar  conclusions  are  also 
reached  in  the  analytical  theories  of  turbulence  (see  the 
monographs  of  Lesieur1  and  McComb18),  which  predict 
that  the  energy  flux  through  a  wave  number  k  in  the  iner¬ 
tial  range  is  dominated  by  scales  from  the  spectral  vicinity 
of  k.  However,  these  conclusions  are  not  necessarily  incon¬ 
sistent  with  the  nonlocal  character  of  the  detailed  transfer 
function  T(k\p,q)  found  by  Domaradzki  and  Rogallo2,3 
since  the  energy  flux  across  the  energy  spectrum  is  ob¬ 
tained  as  a  wave  number  integral  of  T{k\p,q),  and  because 
of  the  cancellation  effects,  the  integral  may  not  reflect  the 
nonlocal  character  of  the  integrand.  One  must  therefore  be 
cautious  in  drawing  conclusions  about  the  detailed  transfer 
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T\k  p.q)  from  the  analyses  of  integrated  quantities. 

Some  ot  these  theories  may  provide  a  partial  explana¬ 
tion  of  features  of  the  transfer  process  observed  in  the  sim¬ 
ulations,  but  none  of  them  is  fully  satisfactory.  As  already 
noted,  the  classical  theories  based  on  the  spectral  energy 
flux  assumptions  do  not  account  explicitly  for  the  nonlocal 
interactions.  On  the  other  hand,  the  linear  theories  of  the 
far  dissipation  range,  which  have  nonlocal  character,  lead 
to  the  form  of  the  energy  spectrum  proportional  to 
exp(  —air),  contrary  to  the  form  exp(  —ak ),  which  is  ob¬ 
served  in  the  dissipation  range  of  low  Reynolds  number 
experiments  (Ling  and  Huang19  and  Comte-Bellot  and 
Corrsin20),  simulations  (Kerr21  and  Kida  and  Mura¬ 
kami"),  and  is  also  predicted  theoretically  by  Foias  et  aL 23 
and  the  analytical  theories  of  turbulence  (Kraichnan24). 
Among  linear  theories  only  Kraichnan’s  analysis15  shows  a 
possibility  of  the  behavior  exp(—  ak)  for  the  spectrum  of 
the  scalar  field  in  the  viscous-conductive  range.  However, 
this  result  is  obtained  under  assumptions  of  unknown  va¬ 
lidity,  and  no  parallel  result  for  the  far  dissipation  range 
behavior  exists.  The  asymptotic  analysis  of  Brasseur8  re¬ 
sults  in  a  scaling  of  the  energy  transfer  rate  between  two 
small  scales  by  the  square  root  of  the  energy  of  the  remain¬ 
ing  large  scale  rather  than  by  its  energy,  as  suggested  by 
the  phenomenological  relation  (7).  Also,  Waleffe’s 
analysis9  is  unable  to  produce  this  scaling.  The  EDQNM 
theory  reproduces  properly  the  scaling  (7)  for  power  law 
spectra  in  the  disparate-scale  limit  (Domaradzki  and 
Rogallo3),  in  particular,  the  dependence  of  T{k\p,q)  on 
the  energy  E{q )  of  the  large  scales.  However,  the  physical 
interpretation  of  the  observed  transfer  mechanism  is  not 
made  easier  by  referring  to  the  EDQNM  theory,  which  is 
derived  using  assumptions  about  statistical  properties  of 
turbulent  fields  and  formal  structure  of  the  Navier-Stokes 
equations  rather  than  an  intuitive  picture  of  the  physical 
processes  occurring  in  turbulent  flows. 

HI,  PHYSICAL  INTERPRETATION  OF  THE  SELF- 
SIMILARITY  SCALING 

Dimensional  analysis  of  the  phenomenological  relation 
(7)  indicates  that  Ts  has  the  dimension  of  time.  Since  the 
function  Ts  is  the  same  for  all  pairs  of  wave  numbers  p  and 
q  it  must  scale  with  a  time  scale  of  turbulence  that  does  not 
depend  explicitly  on  p  and  q.  It  is  thus  natural  to  assume 
that  it  is  a  certain  integral  time  scale  of  turbulence,  e.g., 
large  eddy  turnover  time, 

/  f*c  \-l/2 

©~y  I^Ei^dkj  ,  (8) 

where  kc  denotes  the  wave  number  at  the  end  of  the 
energy-containing  range.  The  quantity  (8)  has  been  com¬ 
puted  for  all  three  velocity  fields  considered  in  this  paper 
and  is  compared  in  Table  1  with  the  values  of  the  positive 
peaks  of  the  function  Ts  for  these  fields.  The  peak  values 
are  within  a  factor  2  of  the  computed  integral  time  scales. 
This  is  probably  as  good  an  agreement  as  could  be  ex¬ 
pected,  in  view  of  the  fact  that  the  formula  (8)  provides 
only  an  order-of-magnitude  estimate  of  the  large  eddy 
turnover  time. 
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TABLE  i  Large  edd>  turnover  times  @  and  the  peak  values  of  the 
self-simiiar  mncnor.  T  for  three  different,  numencal  veiocitv  fields. 


Field 

© 

T  max 

J  i 

K.128 

0.11 

0.230 

C128U8 

0.02 

0.026 

LES128 

0.04 

0.077 

Consider  now  the  energy  transfer  rate  to  all  wave  num¬ 
bers  k>p  interacting  with  the  wave-number  bands  (p,p 
+  A p)  and  (q,q+Aq).  This  transfer  rate  is  equal  to  an  area 
under  the  positive  peak  of  the  function  T(k\p, q)  and  its 
estimate,  using  (7),  is 

TI(p,q)~  X  T(k\p,q)Ap  Aq~qpE(p)E(q)Ap  Aq  9. 

k>p 

(9) 

The  factor  q  in  the  last  equality  appears  because  the 
range  of  wave  numbers  k  >p  affected  by  the  wave-number 
bands p  and  q  is p<,k<,p+q,  as  discussed  after  formula  (7). 
Due  to  the  localness  of  the  energy  transfer  process,  the 
quantity  II  (p,q)  may  also  be  interpreted  as  the  energy 
transfer  rate  from  the  scales  in  the  band  (p,/H-Ap)  caused 
by  their  interactions  with  scales  from  the  bands  (q,q+&q) 
and  ( p,P+q )• 

The  classical  spectral  energy  transfer  hypotheses  usu¬ 
ally  invoke  dimensional  arguments,  eddy  viscosity  con¬ 
cepts,  or  analogies  with  the  turbulent  energy  production  by 
the  mean  flow  gradients  in  the  turbulent  kinetic  energy 
balance  equation.  The  last  approach,  used  by  Obukhov25 
and  Ellison,26  is  the  most  appealing  because  it  is  based  on 
the  formally  correct  equation.  The  main  difficulty  in  this 
approach  is  a  proper  interpretation  of  the  equation  in  terms 
of  spectral  quantities. 

The  turbulent  energy  production  term  is 

- dU, 

P=-U'U^'  (10) 

where  up j  is  the  Reynolds  stress  tensor  and  U,  is  the  mean 
flow.  We  will  attempt  to  express  (9)  in  a  form  consistent 
with  (10)  so  that  the  Reynolds  stresses  and  the  rates  of 
strain  of  the  large-scale  flow  can  be  identified. 

The  following  notation  is  introduced:  Lq=\/q — a 
length  scale  of  large  eddies  from  the  q  band;  lp=\/p — a 
length  scale  of  small  eddies  from  the  p  band; 
U)j=E(q)Ci.q— energy  of  the  large  scales;  up 
=£(p)A/> — energy  of  the  small  scales.  Using  this  notation 
(9)  can  be  rearranged  as  follows: 

T[{p,q)~\E{p)Lp){WE{q)bq)m[p1E{q)Lq\'n<S>} 


The  formula  (11)  has  the  same  form  as  (10)  if  the 
following  identification  is  made: 

(12) 


dU  U,  U. 

- - -py©. 

ox,  Lq  lp 

The  estimate  (12)  of  the  Reynolds  stress  is  consistent 
with  the  classical  arguments  (see  Mornn  and  Yaglom,'  pp. 
212-225),  which,  on  the  basis  of  the  dimensional  consid¬ 
erations,  postulate  proportionality  between  the  Reynolds 
stresses  and  energy.  Departure  from  the  classical  argu¬ 
ments  occurs  in  the  estimate  of  dU/dxr  The  estimate  ( 13 ) 
is  a  composite  expression  that  involve*  a  product  of  the 
classical  estimate  U/Lg  time  scale  of  turbulence  ©,  and 
the  rate  of  strain  Ug/Ig.  Since  V \/Iq>  k’/Lg  and  ©> 
this  product  is  much  larger  than  the  classical  estimate 
Uq/Lg,  and  it  suggests  the  existence  of  flow  regions  with 
the  rates  of  strain  determined  by  the  velocity  scale  of  the 
large  eddies  and  the  length  scale  of  the  small  eddies.  Such 
large  strains  may  be  induced  by  a  mutual  interaction  of  a 
few  large  scales  occurring  over  long  enough  time.  For  in¬ 
stance,  in  boundary  layer  flows  counter-rotating  stream- 
wise  vortices  with  a  length  scale  comparable  to  the 
boundary-layer  thickness  8  are  known  to  generate,  in  about 
one  large  eddy  turnover  time,  internal  shear  layers  with 
much  smaller  length  scales  l <8  and  the  velocity  scale  com¬ 
parable  to  the  mean  ffee-stream  velocity,  i.e.,  the  velocity 
scale  of  the  largest  eddies.  Secondary  instabilities  of  these 
shear  layers  are  responsible  for  the  transition  to  turbulence 
and  the  generation  of  small  scales,  as  demonstrated  exper¬ 
imentally  by  Swearingen  and  Blackwelder27  and  numeri¬ 
cally  by  Liu  and  Domaradzki.28  In  free  shear  layer  flows 
similar  strong  internal  shears  have  also  been  observed. 
Thus  the  existence  of  intermittent  regions  with  strains 
much  larger  than  postulated  by  the  classical  theories  of 
turbulence  is  plausible,  even  though  it  has  not  been  sys¬ 
tematically  investigated  for  isotropic  turbulence.  It  is  in¬ 
teresting  to  note  that  in  the  context  of  a  passive  scalar  in 
the  viscous-conductive  range  Kraichnan15  observes  that 
large  intermittent  rates  of  strain,  if  they  exist  in  a  flow,  will 
determine  the  dissipation  spectrum  of  the  scalar.  The 
above  analysis  of  the  scaling  (7)  provides  a  plausible  phys¬ 
ical  mechanism  of  the  energy  transfer  process,  where  the 
role  of  the  large  scales  is  to  produce  intermittent  regions  of 
relatively  strong,  internal  shears  characterized  by  smaller 
length  scales,  which  serve  as  regions  of  efficient  small-scale 
energy  transfer.  If  one  notes  that  ©  represents  the  inverse 
of  a  collective  rate  of  strain  of  all  scales  from  the  energy- 
containing  range  then  (13)  is  a  fraction  of  the  total  strain 
Ug/lg  which  is  attributable  to  the  action  of  the  large  scales 
from  the  wave-number  band  {q,q-±-Lq).  The  development 
of  such  intermittent  regions  is  the  result  of  an  evolution  of 
a  flow  over  about  one  large  eddy  turnover  time,  and  it  is 
unlikely  that  any  analysis  that  does  not  explicitly  account 
for  time  evolution  will  be  able  to  predict  scaling  (7).  Note 
also  that  in  this  interpretation  dependence  of  the  transfer 
on  rates  of  strain  in  the  physical  space  is  recovered,  despite 
its  formal  dependence  on  the  large-scale  energy  in  the  spec¬ 
tral  representation. 

To  confirm  (or  reject)  the  interpretation  of  the  energy 
transfer  process  proposed  above  a  detailed  analysis  of  tur¬ 
bulent  fields  should  be  made. 
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IV.  CONSEQUENCES  OF  THE  SCALING  RELATION 
FOR  LOW  REYNOLDS  NUMBER  FLOWS 


Since  the  classical  closure  hypotheses  usually  concern 
the  energy  flux  through  a  given  wave  number,  say  k. 


W(k)  = 


T{k')dk'  = 


(14) 


i.e.,  the  energy  transfer  rate  o  all  wave  numbers  k'>k 
from  all  wave  numbers  k'  <,k,  it  is  important  to  investigate 
the  consequences  of  the  scaling  (11)  on  the  flux  (14). 

The  flux  (14)  is  obtained  from  (11)  by  considering 
transfer  from  all  p<,k  and  q<,k.  Because  of  the  locality  of 
the  energy  transfer  only  wave  numbers  p  in  the  vicinity  of 
k,  for  instance  between  1/2 k  and  k,  will  contribute.  There¬ 
fore. 


fV(k)~  f  dp  f  dqU(p,q) 

J  ( 1/2)*  Jo 

~0(/(  ^pE(p)dp^  qE(q)dq^.  (15) 


au 

dx, 


!  r*  ■  ’-sz  it 

k'-El  k'  idk'  j  ~  — 

■Jo  L 


The  closure  hypothesis  (20)  may  be  used  to  predict  the 
form  of  the  energy  spectrum  for  the  wave  numbers  beyond 
the  energy-containing  range,  where  it  applies.  The  spectral 
energy  equation  is 


jt  j*E(k',t)dk'  =  -JV(k,t)-2v ^  k'2E(k')dk'. 

(22) 

In  turbulence  at  low  Reynolds  numbers  the  energy 
spectrum  falls  off  rapidly  for  increasing  wave  numbers  k 
and  for  the  wave  numbers  beyond  the  energy-containing 
range, 


j*  E(k’,t)dk'~  j* 


E(k’,t)dk'  u2(t),  (23) 


where  u  is  the  root  mean  square  of  turbulent  velocity. 

Using  the  velocity  scale  U(k)—u,  Eq.  (22)  may  be 
rewritten  as 


For  k  beyond  the  energy-containing  range  the  integrals 
in  (15)  can  be  estimated  as  follows: 


€=$frE(k)u-t-2v 


k'lE{k’)dk' , 


(24) 


f  pE(p)dp~k2E(k), 

1  (1/2)* 

(16) 

rk  U 2  U 

J  oqE(q)dq~.T~-, 

(17) 

where  U  and  L  are  the  integral  velocity  and  length  scale, 
respectively.  Note  that,  in  general,  these  quantities  depend 
on  the  wave  number  k,  for  instance  U—U{k).  Using  these 
formulas  the  estimate  of  the  energy  flux  becomes 


where  e  is  the  energy  dissipation  rate,  P  is  a  constant,  and 
the  explicit  time  dependence  is  omitted.  Differentiating 
(24)  with  respect  to  k  we  obtain  the  equation 

d 

Pu—[k2E(k)]+2vk1E(k)=0,  (25) 


which  has  the  following  solution: 
E(k)~-pexp(—ak), 


(26) 


W(k)~k?E(k)U{k).  (18) 

The  same  result  is  obtained  more  directly  by  assuming 
that  in  (11)  Uq  and  Lq  are  replaced  by  the  integral  scales 
U  and  L,  and  the  rate  of  strain  of  small  scales  u2  is  ap¬ 
proximated  by  Ellison’s  hypothesis:13 

u2~pE(p).  (19) 

Note,  however,  that  the  derivation  of  Eq.  (18)  relies  only 
on  the  scaling  relation  (7)  and  not  on  its  physical  inter¬ 
pretation  given  in  the  previous  section.  Therefore,  the  re¬ 
sults  of  this  and  the  next  section  are  independent  of  the 
particular  physical  interpretation  proposed  in  this  work. 

In  summary,  the  formula  (18)  corresponds  to  the  fol¬ 
lowing  closure  assumption  for  the  spectral  energy  flux: 

W{k)~-w.d-!^~  [kE(k)  ][*£f(*)],  (20) 

dxj 

where  the  Reynolds  stress  of  the  small  scales  is  represented 
by  (19)  and  the  rate  of  strain  acting  on  these  scales  is  given 
as  a  ratio  of  the  integral  velocity  scale  U  and  the  length 
scale  of  the  small  eddies  l/k.  Note  that  in  both  Obukhov 
and  Ellison  theories  this  rate  of  strain  is  estimated  by  the 
following  expression: 


where  o=  (2v)/(fJu).  Equation  (25)  expresses  an  approx¬ 
imate  balance  between  viscous  dissipation  and  nonlinear 
transfer 

dW 

m)  =  -— ~exp(-a/:)  (27) 

in  the  far  dissipation  range.  Equations  (26)  and  (27)  are 
important  conclusions  derived  from  the  scaling  relation 
(7),  which  can  be  compared  with  experimental  and  nu¬ 
merical  results  for  low  Reynolds  number  turbulence. 

In  Fig.  2(a)  we  plot  the  energy  spectrum  for  isotropic 
field  K128  and  in  Fig.  2(b)  the  spectral  energy  balance  for 
this  field.  At  this  low  Reynolds  number,  25,  all  spec¬ 
tra  are  well  resolved,  with  both  the  dissipation  spectrum 
and  the  transfer  spectrum  vanishing  for  large  wave  num¬ 
bers  k  and  in  the  approximate  balance  outside  the  energy- 
containing  range  0<k<  10,  where  the  above  derived  ex¬ 
pressions  should  be  valid.  Using  the  log-linear  scales  we 
plot  in  Fig.  3  the  dissipation  spectrum,  i.e.,  the  energy 
spectrum  multiplied  by  k2,  and  the  transfer  spectrum  for 
this  field  with  wave  numbers  k  rescaled  using  the  Kolmog¬ 
orov  length  7]  zb  0.0344.  It  is  seen  that  the  functional  fonm. 
Eqs.  (26)  and  (27)  are  in  an  excellent  agreement  with  the 
numerical  results  for  the  wave  numbers  k>10  ( t^A>0.3 ), 
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FIO.  2.  (*)  Energy  spectrum  for  isotropic  field  K128  plotted  using  linear 
scales  to  accentuate  the  energy-containing  range  0<* <  10.  (b)  Spectral 
energy  balance  for  the  field  K128.  Dashed  line,  dissipation  spectrum 
.  —2 vk*E(k).  dash-dotted  line,  nonlinear  energy  transfer  T(k);  dotted 
line,  sum  of  dissipation  and  transfer,  T(k)  —  2vk*E{k)  =dE/dt. 
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FIG.  4.  Dissipation  spectrum  2vk1E(k)  and  transfer  spectrum  T(k) 
beyond  the  energy-containing  range  for  the  turbulent  Taylor-Green  vor¬ 
tex  field  plotted  using  log-linear  scales  and  wave  number  k  normalized  by 
the  Kolmogorov  length  seal'  -j. 

i.e.,  for  all  wave  numbers  beyond  the  energy-containing 
range.  The  slope  of  the  curves  in  Fig.  3  is  a/17  ~  2. 19. 

In  Fig.  4  we  plot  in  the  same  manner  the  energy  and 
transfer  spectra  for  a  flow  at  higher  Reynolds  number, 
Rx~ 70,  obtained  in  numerical  simulations  of  a  decaying 
Taylor-Green  vortex  performed  with  an  effective  resolu¬ 
tion  5123  modes  (Brachet29).  For  the  wave  numbers  be¬ 
yond  the  energy  containing  range  Eqs.  (26)  and  (27)  pro¬ 
vide  a  very  good  fit  with  the  nondimensional  constant 
fl/i?=;2.08. 

Finally,  in  Fig.  5  the  experimental  results  of  Comte 
Bellot  and  Corrsin20  for  the  dissipation  spectrum  of  grid 
turbulence  at  R^  =  60.7  and  corresponding  numerical  re¬ 
sults  of  Ruetsch  and  Maxey30  for  forced  turbulence  at 
•K,.  —  60  are  shown,  again  exhibiting  good  agreement  (ex¬ 
perimental  and  numerical  data  form  a  straight  line  on  the 


FIG.  3.  Dissipation  2vk‘E(k)  (solid  line)  and  transfer  spectrum  T(k) 
beyond  the  energy-containing  range  (broken  line)  for  the  field  K128 
plotted  usmg  log-linear  scales  and  wave  number  k  normalized  by  the 
Kolmogorov  length  scale  ij. 


FIG.  5  The  normalized  dissipation  spectrum  2vtcE(k)  in  experiments  of 
Comte-Bellot  and  Cornin'*  (O)  and  numerical  simulations  of  forced 
isotropic  turbulence  of  Ruetsch  and  Maxey1*  (  x  ). 
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log-linear  plot'  with  formula  >26)  beyond  the  energy- 
containing  range.  The  nonaimensionai  consiam  a; -77  — 2.12 
in  this  plot. 

In  view  of  the  qualitative  (scaling)  character  of  the 
assumptions  leading  to  the  closure  (18)  the  agreement  be¬ 
tween  predictions  ( 26 )  and  ( 27 )  and  the  observed  form  of 
the  energy  and  transfer  spectra  is  surprisingly  good.  In  a 
similar  range  of  Reynolds  numbers  Kerr1  fitted  numerical 
energy  spectra  by  a  formula  E{k)~k~in  exp(  —  ak),  with 
the  nondimensional  constant  <2/17  —  5.1.  and  Kida  and 
Murakami2*  obtained  the  best  fit  to  the  energy  spectrum  in 
the  dissipation  range  as  k~]  0  exp(  —  ak)  with  a/p~ 4.9. 
Sreenivasan31  used  as  the  best  fit  to  the  experimental  en¬ 
ergy  spectra  the  formula  exp t -ak)  with  <2/77=  12.7  for 
0.1<ki7<0.5  and  <2/17  =  8.8  for  Q.5<kq<,\.5.  These  results 
and  Eq.  (26)  all  predict  the  existence  of  the  exponential 
factor  in  the  formula  for  the  dissipation  range  of  low  Rey¬ 
nolds  number  turbulence,  but  they  differ  in  the  form  of  the 
algebraic  prefactor.  At  the  present  time  the  quality  of  the 
available  data  is  probably  not  sufficient  to  distinguish  be¬ 
tween  different  prefactors.  In  any  case,  the  exact  value  of 
the.  exponent  m  in  the  prefactor  proportional  to  (krj)~m 
with  positive  m  will  influence  the  behavior  of  the  energy 
spectrum  only  for  kp=0(\).  Note,  however,  that  the 
value  of  the  exponent  <2/77  will  depend  on  the  exact  form  of 
the  prefactor,  decreasing  for  increasing  m,  since  for  a  given 
energy  spectrum  the  exponential  function  must  compen¬ 
sate  for  different  fall-off  rates  associated  with  different  pre¬ 
factors. 

At  this  point  it  is  appropriate  to  briefly  compare  our 
approach  to  compute  the  energy  flux  with  the  recent  re¬ 
sults  of  Waleffe.9  His  analysis,  performed  for  the  inertial 
range  spectrum,  leads  to  the  conclusion  that  two  distinct 
classes  of  mad  interactions  operate  in  turbulent  flows.  Tri¬ 
ads  in  class  R  (for  reverse  transfer)  are  characterized  by 
the  middle  leg  in  the  triad  transferring  large  amounts  of 
energy  to  the  longest  leg  (small  scale)  and  small  amounts 
to  the  shortest  leg  (large  scale).  In  triads  belonging  to  the 
class  F  (for  forward  transfer)  the  shortest  wave  number  is 
losing  energy  to  the  middle  and  the  smallest  wave  num¬ 
bers.  The  triads  considered  as  dominant  in  the  present 
work,  which  are  nonlocal  but  result  in  the  local  energy 
transfer,  belong  to  class  R  in  Waleffe’s  nomenclature.  An 
important  result  from  his  analysis  is  that  the  energy  flux 
(14)  is  determined  primarily  by  class  F  interactions  and 
the  weaker  effect  of  class  R  interactions  is  the  inverse  en¬ 
ergy  flux,  from  small  to  large  scales. 

In  view  of  this  result  our  approach  to  compute  the 
energy  flux  using  Eq.  (15),  which  takes  into  account  only 
the  class  R  triads,  should  be  incorrect.  To  resolve  this 
apparent  contradiction  between  the  present  approach  and 
Waleffe’s9  theory  we  have  calculated  for  the  isotropic  field 
K128  the  energy  flux  using  Eq.  (14)  and  the  approxima¬ 
tion  analogous  to  (15),  where  the  numerically  computed 
T{k'\p,q)  is  integrated  over  0 <q<\k,  \k<p<k,  and 
k'  >  k.  The  inequalities  for  p  and  q  also  imply  that  the 
length  of  the  third  leg  of  the  triad  does  not  exceed  }k.  Both 
results  are  compared  in  Fig.  6  for  wave  numbers  k  >  20 
(17/1=0.688),  i.e.,  outside  the  energy-containing  range. 
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FIG.  6.  Total  energy  riux  k ,  isoiid  line)  and  contribution  to  the  flux 
from  tnads  with  wave  numbers  0 < q < t/c.  < p  <  k,  and  k<k’  <\k  ( bro¬ 
ken  line).  The  axes  are  scaled  as  in  Fig.  3  to  facilitate  comparison  of 
slopes  in  both  figures. 

Clearly,  the  approximation  is  in  an  excellent  agreement 
with  the  exact  energy  flux,  and  we  thus  conclude  that  it  is 
fully  supported  by  the  numerical  data  in  the  dissipation 
range.  Moreover,  formula  (18)  and  the  predicted  form  of 
the  energy  spectrum  (26)  imply  that  the  energy  flux 
should  be  proportional  to  exp  (—ak).  The  numerical  re¬ 
sults  plotted  in  Fig.  6  are  in  agreement  with  this  conclu¬ 
sion. 

It  may  also  be  noted  that  as  seen  in  Fig.  6  these  R-type 
interactions  generally  overestimate  the  total  energy  flux 
(but  by  no  more  than  10%).  We  may  thus  expect  that  the 
remaining  triads,  mostly  of  the  F  type,  will  result  in  the 
inverse  energy  transfer  to  counter  the  excess  of  the  energy 
transferred  to  small  scales  by  the  R  interactions.  Therefore, 
the  numerical  results  point  to  the  energy  transfer  mecha¬ 
nism  in  the  dissipation  range,  large  forward  transfer  by  the 
class  R  tnads,  and  small  inverse  transfer  by  the  F  class 
interactions,  which  is  in  the  diametrical  opposition  to  the 
mechanism  proposed  by  Waleffe9  for  the  inertial  range  dy¬ 
namics.  Among  other  numerical  works,  which  also  do  not 
agree  with  Waleffe’s9  theory  are  the  results  of  Ohkitani  and 
Kida.6  Those  authors  found  that  in  their  simulations  of  a 
flow  at  R^~  1 80,  three  types  of  interactions  dominated. 
Among  them,  the  class  R  triads  constituted  about  30%  of 
all  active  triads,  the  class  F  about  15%,  and  more  than 
50%  triads  were  characterized  by  both  the  largest  and  the 
middle  scale  losing  energy  to  the  smallest  one,  thus  not 
adhering  to  Wallefe’s  classification.  In  view  of  the  above 
disagreements  it  seems  evident  that  the  turbulence  dynam¬ 
ics  at  these  lower  Reynolds  numbers  m3'  >  substantially 
different  from  the  dynamics  of  turbuler  r  the  inertial 
range. 

V.  CONSEQUENCES  OF  THE  SCALING  RELATION 
FOR  HIGH  REYNOLDS  NUMBER  FLOWS 

Closure  (20)  is  based  on  the  phenomenological  rela¬ 
tion  (7),  which  was  obtained  from  the  results  of  the  direct 
numerical  simulations  of  low  Reynolds  number  turbulence. 
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However,  relation  <")  was  also  found  to  hold  for  the  ve¬ 
locity  field  obtained  in  a  large  eddy  simulation  of  turbu¬ 
lence  at  nominally  infinite  Reynolds  number  (LES12S).  It 
may  therefore  be  useful  to  formally  explore  the  conse¬ 
quences  of  this  closure  for  high  Reynolds  number  flows. 

In  the  far  dissipation  range  of  high  Reynolds  number 
turbulence  Eq.  (22)  still  holds,  but  it  is  unlikely  that  the 
velocity  scale  will  be  given  by  the  formula  (23),  which 
emphasizes  scales  from  the  energy-containing  and  inertial 
ranges.  In  this  case  we  follow  Batchelor,14  Novikov,16  and 
Saffmann,17  who  assumed  that  the  far  dissipation  range 
eddies  interact  most  strongly  with  the  Kolmogorov  eddies 
characterized  by  the  velocity  scale  vk=  (ve)i/4.  With  this 
velocity  scale  the  formula  (26)  becomes 


£U)~pexp^  -jj  (r ]k)  ), 


(28) 


where  r;  is  the  Kolmogorov  length  scale: 

The  exponential  factor  in  (28)  is  in  agreement  with 
some  other  theories  of  the  far  dissipation  range.  For  in¬ 
stance,  Kraichnan  predicts  E(k)  —/r3  exp ( —ak)  in  the 
framework  of  the  direct  interaction  approximation  theory 
and  Foias  et  aL2i  obtain  an  exact  sharp  estimate  E(k) 
=o[exp(  —  crjk)]  using  spectral  properties  of  the  Stokes 
operator.  That  work  also  suggests  the  prefactor  k~ 4,  but 
this  result  is  not  considered  mathematically  exact  by  the 
authors.  The  faster  falloff,  proportional  to  exp(  —bk2),  is 
suggested  by  a  number  of  older  theories  (see  Monin  and 
Yaglomb)  and  the  more  recent  work  of  Smith  and 
Reynolds.32  However,  Manley33  showed  that  the  procedure 
employed  by  Smith  and  Reynolds,32  after  a  reasonable 
modification,  provides  results  in  agreement  with  the 
exp  {—ak)  behavior.  It  seems  that,  at  the  present  time,  the 
experimental  and  numerical  results  favor  the  exp  {—ak) 
behavior  with  a  simple,  algebraic  prefactor,  whose  exact 
form  remains  to  be  determined. 

In  the  inertial  range  the  energy  spectrum  decays  too 
slowly  to  use  the  approximation  (23).  We  take  instead 


tf-Cf(*,f)~(  ^E(k',t)dk' 


1/2 


(29) 


To  avoid  divergence  of  the  integral  at  small  k  the  in¬ 
ertial  form  of  the  energy  spectrum  cannot  be  continued  to 
k=0  and  the  specific  form  of  the  energy  spectrum,  which 
vanishes  at  the  low  wave  numbers,  must  be  used. 

Because  of  the  limitation  to  the  inertial  range  wave 
numbers  the  viscous  term  in  Eq.  (22)  may  be  neglected, 
and  using  (29)  one  obtains  the  following  equation: 


f=0krE(k) 


f  *  E{k',t)dk' 
Jo 


1/2 


(30) 


Introducing  a  variable  8(k)  =S^E(k' ,t)dk' ,  Eq.  (30) 
is  reduced  to 


,  ,„dd 
e=pk~d  dk' 

The  solution  to  (31)  is 


(31) 


0U;=- u:,  1-j  '.32) 

where  u  =  []S*E(k',t)dk']l/2  is  the  mean  turbulent  veloc¬ 
ity.  For  decaying  turbulence  e  —  u}/L,  where  L  is  the  inte¬ 
gral  length  scale  of  turbulence.  Then,  if  the  constant  0  is  of 
the  order  unity  and  k>\/L  (inertial  range),  we  may  ex¬ 
pand  the  expression  (32)  in  Taylor  series  and  keep  only 
the  lowest-order  term  in  1  /k.  Differentiating  this  result 
with  respect  to  k  we  obtain  the  energy  spectrum  for  the 
wave  numbers  1/Z-<A:< I/77: 


E(k)  = 


1/3 


e273*- 


5/3 


(33) 


i.e.,  the  classical  Kolmogorov  inertial  range  form.  Assum¬ 
ing  the  value  of  the  Kolmogorov  constant  CK=r  1. 5-2.0,  we 
obtain  for  the  constant  0  a  range  of  values  0.44-0.29. 

The  coefficient  a/r)=2/0  in  the  exponential  law  (28) 
for  this  range  of  values  of  the  Kolmogorov  constant  is 
between  4.5  and  6.9.  This  value  does  not  agree  with  the  low 
Reynolds  number  result  a/17 =2.  The  disagreement  may 
reflect  differences  between  low  and  high  Reynolds  number 
flows  or,  what  is  more  plausible,  an  inadequacy  of  the  as¬ 
sumption  that  the  far  dissipation  range  eddies  interact 
most  strongly  with  the  eddies  characterized  by  the  velocity 
scale  u=vk=(ve)w*.  Indeed, 


a  2v  1  2  vk 

17  0uij  0  u  ’ 


(34) 


and  if  the  velocity  scale  u  is  replaced  by  the  physically 
more  relevant  velocity  scale  associated  with  the  peak  of  the 
dissipation  spectrum,  vd>vk,  the  constant  a/17  wih  de¬ 
crease.  However,  it  does  not  seem  that  the  value  of  the 
constant  a/17  for  the  far  dissipation  range  can  be  uniquely 
determined  in  the  framework  of  the  above  scaling  argu¬ 
ments. 


VI.  CONCLUSIONS 

Employing  the  phenomenological  scaling  relation  for 
the  detailed  energy  transfer  functions  T(k\p,q)  the  physi¬ 
cal  interpretation  of  the  observed  local  energy  transfer 
caused  by  nonlocal  triads  is  proposed.  In  this  interpreta¬ 
tion  the  large  scales  of  turbulence,  through  their  mutual 
interactions,  create  intermittent  regions  of  relatively  large 
rates  of  strain,  where  most  of  the  small-scale  energy  trans¬ 
fer  occurs.  An  important  feature  of  this  interpretation,  sug¬ 
gested  by  the  presence  of  the  large  eddy  turnover  time  in 
the  scaling  relation,  is  an  introduction  of  the  time  scale  into 
the  problem.  This  implies  that  any  successful  explanation 
of  the  observed  energy  transfer  process  must  invoke  the 
time  evolution  of  turbulence.  It  may  be  an  important  ob¬ 
servation  in  view  of  attempts  to  explain  the  transfer  process 
by  analyzing  Navier-Stokes  equations  considered  at  a 
given  instant  of  time.  In  such  analyses  valid  assumptions 
about  amplitudes  of  Fourier  modes  are  made  and  conclu¬ 
sions  about  the  dependence  of  energy  transfer  on  the  am¬ 
plitudes  are  drawn.  The  transfer  process,  however,  is  de¬ 
termined  by  both  the  amplitudes  and  the  phases  of  the 
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FGurier  modes  and  the  phase  relationships  are  not  known 
a  priori  and  must  be  found  by  considering  the  time  evolu¬ 
tion  of  the  system.  This  fact  may  explain  why  the  scaling 
relation  (7)  is  broadly  consistent  with  the  analytical  theo¬ 
ries  of  turbulence,  which  take  into  account  the  time  evo¬ 
lution  of  turbulence,  and  why  other  approaches  failed  to 
predict  this  scaling.  In  this  respect  it  must  be  stressed  that 
the  scaling  relation  (7)  may  be  treated  as  being  equivalent 
to  an  experimental  result,  and  thus  any  theoretical  expla¬ 
nation  of  the  energy  transfer  process  should  be  required  to 
predict  it. 

Independently  of  its  physical  interpretation  the  scaling 
relation  (7)  was  used  to  deduce  the  functional  forms  of  the 
energy  and  transfer  spectra  in  the  far  dissipation  range  of 
low  Reynolds  number  turbulence.  Considering  an  approx¬ 
imate  character  of  the  analysis,  the  predicted  forms  are  in 
surprisingly  good  agreement  with  several  numerical  and 
experimental  results  for  the  range  of  nondimensional  wave 
numbers  krj  <  2. 

Assuming  the  validity  of  the  scaling  relation  for  high 
Reynolds  number  flows,  the  k~ 5/3  behavior  of  the  energy 
spectrum  is  predicted  in  the  inertial  range  and  the 
k~2  exp(  —ak)  form  in  the  far  dissipation  range.  These 
results  are  on  a  less  Arm  ground  than  the  results  obtained 
for  low  Reynolds  number  flows.  The  reason  is  that  the 
applicability  of  the  scaling  relation  is  confirmed  only  for 
large  eddy  simulation  results  of  high  Reynolds  number 
turbulence,  which  are  influenced  by  the  additional  assump¬ 
tions  made  in  this  numerical  methodology.  Therefore,  the 
more  appropriate  conclusion  would  be  that  the  scaling  re¬ 
lation  (7)  does  not  contradict  the  inertial  range  and  the  far 
dissipation  range  behavior  of  high  Reynolds  number  tur¬ 
bulence. 
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Abstract 


Using  a  velocity  field  obtained  in  a  direct  numerical  simulation  of  isotropic  turbu¬ 
lence  at  a  moderate  Reynolds  number  we  analyze  the  subgrid-scale  energy  transfer  in 
the  spectral  and  the  physical  space  representation.  The  subgrid-scale  transfer  is  found 
to  be  composed  of  a  forward  and  an  inverse  transfer  components,  both  being  significant 
in  dynamics  of  resolved  scales.  Energy  exchanges  between  the  resolved  and  unresolved 
scales  from  the  vicinity  of  the  cutoff  wave  number  dominate  the  subgrid-scale  processes 
and  the  energetics  of  the  resolved  scales  are  unaffected  by  the  modes  with  wave  numbers 
greater  than  twice  the  cutoff  wave  number.  The  dominance  of  nonlinear  interactions 
among  the  largest  scales  in  the  subgrid-scale  energy  transfer  process  suggests  that  the 
resolved  nonlinear  term  may  serve  as  a  basis  of  a  new  approach  to  the  subgrid-scale 


modeling. 


1  Introduction 


Three  approaches  used  in  numerical  predictions  of  turbulent  flows  are  direct  numerical  simu¬ 
lations  (DNS),  large  eddy  simulations  (LES),  and  Reynolds  averaged  Navier-Stokes  (RANS) 
simulations.  With  currently  available  computer  capabilities  the  applicability  of  the  DNS 
methods  is  limited  to  low  Reynolds  number  turbulence.  In  practical  applications  for  high 
Reynolds  number  flows  the  RANS  techniques  are  used  most  frequently.  The  main  draw¬ 
back  of  this  method  is  the  need  for  introduction  of  a  number  of  phenomenological  closure 
assumptions  and  empirical,  flow  dependent  constants. 

i  o 

The  LES  techniques,  reviewed  by  Rogallo  and  Moin  1  and  more  recently  by  Lesieur  , 
are  a  compromise  approach  between  DNS  and  RANS.  In  the  LES  large,  resolved  scales  of  a 
turbulent  flow  are  simulated  directly,  akin  to  the  DNS  approach,  and  their  interactions  with 
the  small,  unresolved  scales  are  modeled  like  in  the  RANS  approach.  However,  contrary 
to  the  RANS,  only  a  part  of  the  nonlinear  interactions  is  modeled  in  the  LES,  and  since 
the  modeled  interactions  involve  small  scales  (usually  in  the  inertial  range  of  turbulence) 
which  have  more  universal  character  than  flow  dependent  large  scales,  the  hope  is  that  such 
modeling  can  be  accomplished  with  less  empiricism  and  with  greater  help  from  the  theories 
of  homogeneous  turbulence  than  it  is  the  case  for  the  RANS  approach.  At  the  present 
time  the  most  widely  used  subgrid-scale  models  are  the  Smagorinsky  model  °  for  the  LES 
performed  in  the  physical  space  representation  and  the  Kraichnan  ^  and  the  Cholet  and 
Lesieur  models  ^  if  the  spectral  representation  is  used.  These,  as  well  as  other  subgrid-scale 
models,  despite  exhibiting  a  number  of  desirable  properties  like  accounting  properly  for  the 
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global  energy  flux  from  the  large  to  the  small  scales,  are  known  to  be  deficient  in  some 
respects.  For  instance,  the  models  are  usually  purely  dissipative.  However,  the  process  of 
the  subgrid-scale  energy  transfer  is  dissipative  only  in  the  mean  and  locally  in  the  spectral 
or  the  physical  space  the  effect  of  the  subgrid-scale  interactions  may  be  to  either  decrease  or 
increase  the  energy  of  the  large,  resolved  scales.  Various  attempts  were  proposed  in  the  past 
to  account  for  the  inverse  energy  transfer  in  the  subgrid-scale  modeling  for  homogeneous 
turbulence  D,‘’°  and  for  inhomogeneous  flows  y  but  no  generally  accepted  method  exists. 

The  practical  importance  of  the  LES  techniques  and  the  deficiencies  of  the  existing 
subgrid-scale  models  suggest  that  better  understanding  of  subgrid-scale  interactions  is  needed 
if  improvements  in  the  LES  methods  are  to  be  made.  To  compute  the  effects  of  the  subgrid- 
scale  nonlinear  interactions  a  full  velocity  field  in  three-space  dimensions  must  be  known; 
such  detailed  information  cannot  be  obtained  using  current  experimental  techniques.  Re¬ 
quired  information,  however,  is  available  in  the  direct  numerical  simulations  of  turbulent 
flows  and  has  been  used  in  the  past  to  investigate  the  subgrid-scale  interactions  and  to  as¬ 
sess  directly  the  validity  of  the  models.  Such  an  approach  was  pioneered  by  Clark  et  al.  ^ 
for  the  physical  space  modeling  and  by  Domaradzki  et  al.  ^  for  the  spectral  space  mod¬ 
eling.  A  major  limitation  of  this  approach  is  that  only  low  Reynolds  number  flows  can  be 
simulated  numerically  and  thus  it  is  unclear  to  what  extent  conclusions  from  such  analyses 
are  applicable  to  more  important  case  of  high  Reynolds  number  turbulence. 

In  this  work  we  investigate  the  properties  of  the  subgrid-scale  nonlinear  interactions 
using  both  the  physical  and  the  spectral  space  representation  for  numerically  simulated, 
decaying  homogeneous  turbulence.  The  simulated  flow  is  the  Taylor-Green  vortex  and  using 
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i  c\ 

its  symmetries  it  is  possible  to  increase  Reynolds  number  by  a  factor  2  as  compared  with 
the  general  nonsymmetric  flows  simulated  with  the  same  number  of  computational  modes. 
It  is  hoped  that  the  higher  Reynolds  number  and  the  existence  of  a  short  inertial  subrange 
for  this  flow  can  make  results  of  such  an  investigation  applicable  to  high  Reynolds  number 
turbulence. 

2  Numerical  simulations 

The  Taylor-Green  vortex  flow  °  develops  from  the  following  initial  condition: 

u  =  sin(x)cos(y)cos(z) 

v  =  —  cos(x)  sin(y)  cos(z)  (1) 

w  =  0 

At  time  t  =  0  the  flow  is  two-dimensional  but  becomes  three-dimensional  for  all  times  t  >  0 
when  it  develops  into  initially  well  organized,  laminar  structures  in  the  form  of  vortex  sheets 
which  subsequently  become  unstable  resulting  eventually  in  a  fully  turbulent  flow.  It  was 
noted  by  Orszag  that  the  initial  condition  (1)  has  a  number  of  symmetries  which  are 
consistent  with  symmetries  of  the  Navier-Stokes  equations  and  are  thus  preserved  in  time 
as  flow  evolves.  In  the  context  of  spectral  simulations  the  symmetries  of  the  flow  may 
then  be  used  to  reduce  number  of  computational  modes  needed  to  describe  the  flow  for  a 
prescribed  range  of  resolved  scales.  This  idea  was  implemented  by  Brachet  et  al.  ^  who 
were  able  to  simulate  the  Taylor-Green  vortex  flow  with  an  effective  spatial  resolution  of 


25 63  modes  at  a  computer  cost  equivalent  to  simulating  a  general,  non-symmetric  flow  with 
the  resolution  of  643  modes.  More  recently  Brachet  ^  reported  results  of  simulations  of  the 
Taylor-Green  flow  performed  with  an  effective  resolution  of  up  to  8643  modes  and  Reynolds 
number  R\  %  140.  A  similar  approach  to  increase  range  of  scales  and  Reynolds  numbers 
in  numerical  simulations  of  turbulence  by  employing  symmetries  of  Navier-Stokes  equations 
was  pursued  by  Kida  and  his  collaborators  in  a  n-  tber  papers^’^”’^,  employing  a  flow 
with  even  greater  number  of  symmetries  than  the  Taylor-Green  vortex.  At  the  present  time 
these  highly  symmetric  flows  are  the  most  computationally  efficient  means  of  numerically 
simulating  isotropic  turbulence  with  Reynolds  number  R\  on  the  order  100. 

Using  the  numerical  code  developed  by  Brachet  ^  we  have  performed  direct  numerical 
simulations  of  the  Taylor-Green  vortex  flow  in  order  to  generate  a  turbulent  velocity  field 
for  the  purpose  of  an  analysis  of  the  subgrid-scale  nonlinear  interactions.  Since  the  details  of 
such  simulations  were  extensively  described  by  Brachet  et  al.  ^  and  Brachet  ^  we  report 
here  only  a  few  main  features  of  the  time  evolution  of  the  flow  and  its  properties  at  the  end 
of  the  run.  The  velocity  field  at  the  end  of  the  run  is  used  in  the  subsequent  sections  for  the 
analysis  of  the  subgrid-scale  interactions. 

The  flow  is  contained  in  a  cube  with  a  side  length  2x  resulting  in  wave  numbers  k  = 
(ki,  k3)  in  the  spectral  space  with  integer  components  In  the  physical  space  the  flow 
is  periodic  with  the  period  2 ir  in  each  coordinate  direction  x,  y,  and  z.  Because  of  the 
symmetries  the  flow  never  crosses  the  boundaries  x,y,  and  z  =  ir  and  in  the  subsequent 
discussion  it  will  be  visualized  in  the  so-called  impermeable  box  ^  0  <  x,y,z  <  ir.  The 
effective  spatial  resolution  in  the  simulations  was  51 23  modes,  which,  after  dealiasing  by  the 
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2/3  rule,  provides  the  maximum  wave  number  km  =  170  in  each  coordinate  direction.  Since 
the  velocity  and  the  length  scale  of  the  initial  flow  are  order  unity  the  large  eddy  turnover 
time  is  also  order  unity  and  Reynolds  number  is  equal  to  the  inverse  of  molecular  viscosity 
l/r'  (3000  in  the  simulations).  The  simulations  were  run  until  maximum  time  tm  =  18,  i.e. 
for  several  large  eddy  turnover  times,  with  the  time  step  A t  =  0.0025. 

In  Figs.  1(a)  -  1(d)  we  plot  the  time  evolution  of  the  total  turbulent  energy,  the  total 
dissipation  rate  e,  skewness  S,  and  microscale  Reynolds  number  R\,  respectively.  Until  time 
t  «  5  the  evolution  of  the  flow  is  essentially  inviscid  with  the  total  energy  nearly  constant. 
During  this  period  small  scales  are  generated  from  the  initial  condition  (1)  resulting  in  a 
subsequent  rapid  rise  of  the  dissipation  rate  which  peaks  at  t  as  10  and  later  decays  because 
of  the  decrease  in  the  intensity  of  turbulence  caused  by  the  viscous  damping.  The  skewness, 
after  fairly  chaotic  behavior  until  t  as  10  ,  at  the  end  of  the  run  approaches  —0.5,  which  is  the 
generally  accepted  value  for  this  quantity  in  fully  developed  isotropic  turbulence.  The  initial 
value  of  R\  exceeds  1000,  decays  rapidly  becoming  an  order  of  magnitude  less  at  the  time 
of  the  peak  in  the  dissipation  rate  t  as  10,  and  slowly  approaches  the  final  value  R\  ss  70  at 
the  end  of  the  run. 

The  unnormalized  energy  and  dissipation  spectra  at  the  end  of  the  run  are  plotted  m 
Fig.  2(a).  Small  number  of  modes  in  the  low  wave  number  shells  causes  relatively  large 
fluctuations  in  these  quantities  at  low  wave  numbers.  In  the  range  of  wave  numbers  k  <  Jti 
the  energy  spectrum  conforms  to  the  inertial  k~ 5//3  law  with  the  Kolmogoroff  constant  it, 
the  range  2.2  to  2.7.  The  dissipation  spectrum  peaks  at  k  as  20  which,  for  the  calculated 
Kolmogoroff  length  scale  rj  =  0.011  in  the  units  used,  corresponds  to  rjk  «  0.2.  This  value 


agrees  with  experimental  findings  ^  locating  the  dissipation  peak  in  high  Reynolds  number 
turbulence  at  a  wave  number  order  of  magnitude  less  than  1/t?.  Because  of  a  significant 
overlap  of  the  energy  containing  range  and  the  dissipation  range  it  is  unclear  if  the  observed 
inertial  range  spectrum  for  k  <  20  is  the  result  of  the  same  dynamical  processes  that  operate 
at  very  high  Reynolds  numbers  where  there  exists  wide  separation  between  the  energy  and  the 
dissipation  range.  Also  unusually  high  value  of  the  Kolmogoroff  constant  in  the  simulations 
casts  doubt  on  the  significance  of  the  observed  inertial  subrange  as  being  indicative  of  high 
Reynolds  number  turbulence  dynamics.  We  may  claim  at  best  that  the  Reynolds  number 
in  the  simulations  is  high  enough  to  capture  the  beginnings  of  the  inertial  range  dynamics 
but  too  low  to  separate  it  from  the  effects  of  the  dissipation  range  dynamics.  In  the  far 
dissipation  range  for  A:  >  20  the  dissipation  and  the  nonlinear  transfer  spectra  balance  each 
other  and  have  the  functional  form  proportional  to  fc-2exp(— ak)  as  seen  in  Fig.  2(b).  This 
form  was  derived  t :  Uomaradzki  using  scaling  properties  of  the  detailed  energy  transfer 
observed  in  low  Reynolds  number  turbulence. 

3  Basic  quantities 

For  homogeneous  turbulence  incompressible  Navier-Stokes  equations  in  spectral  (Fourier) 
representation  are: 

^un(k)  =  -uk2un(  k)  +  Nn(  k).  (2) 

Here,  ttn(k)  is  tue  velocity  field  in  spectral  space,  with  the  explicit  dependence  on  time 
omitted,  u  is  the  kinematic  viscosity,  and  Afn(k)  is  the  nonlinear  term 
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(3) 


W„(k)  =  -  l-Pnim( k)  J  d3pui(p)um{k  -  p)  , 

where  tensor  Pnfm(k)  accounts  for  the  pressure  and  incompressibility  effects.  The  summation 
convention  is  assumed  throughout. 

Let’s  assume  that  the  wave  number  space  is  divided  into  two  non-overlapping  regions,  £ 
(|k|  <  kc)  signifying  large  scales,  and  S  (|k|  >  kc)  signifying  small  scales.  In  the  LES  ter¬ 
minology  these  scales  are  also  referred  to  as  the  resolved  and  unresolved  scales,  respectively. 
In  the  LES  an  evolution  equation  for  the  velocity  field  un(k)  truncated  to  the  region  £ 


u„(k)  if  ke£ 


(4) 


0  otherwise 

is  sought.  The  truncation  operation  is  trivially  applied  to  the  linear  terms  in  Eq.  (2).  The 
nonlinear  term  (3)  is  decomposed  as  follows.  First,  it  is  computed  with  one  of  the  contributing 
velocity  fields  truncated  to  U  and  the  other  to  V,  where  U  and  V  may  be  any  of  the  two 
previously  prescribed  regions.  Details  of  such  calculations  are  provided  by  Domaradzki  and 
Rogallo  21 ,22  Resuiting  quantity,  denoted  by  Nnuv(k),  describes  the  modification  of  the 
mode  k  caused  by  all  triad  interactions  involving  k  and  two  other  scales,  one  belonging  to  U 
and  the  other  to  V.  Second,  to  retain  the  effect  of  such  nonlinear  interactions  on  the  large 
scales  only,  the  quantity  fVnwv(k)  is  truncated  to  the  region  £,  with  the  result  denoted  by 
NnCUV(k).  The  evolution  equation  for  the  large  scales  £  becomes: 


^UnC(  k)  =  -»k2UnC(  k)  +  ^n(k|*c)  +  Nn’(k\kc), 


(5) 
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where  the  resolved  nonlinear  term  is 

Nn(k\kc)  =  Nnccc(  k),  (6) 

and  the  subgrid-scale  nonlinear  term  Nns{k\kc)  is 

Nna(k\kc)  =  Nnccs(  k)  +  Nncss(k).  (7) 

In  practice,  the  most  straightforward  way  to  compute  (6)  and  (7)  is  to  first  use  (3)  with 
the  full  velocity  fields  u;(p)  and  um(k  —  p)  and  truncate  the  result  to  the  region  £  to  obtain 
the  total  nonlinear  term 


Nntot(k\kc)  =  Nn(k\kc)  +  Nn'(k\kc)  (8) 

Next,  Eq.  (3)  is  used  again  with  the  truncated  velocity  fields  uf( p)  and  u£(k  —  p)  and 
the  result  is  truncated  to  the  region  £  giving  the  resolved  nonlinear  term  ;Vn(k|fcc)  (Eq. 
(6)).  The  subgrid-scale  nonlinear  term  (7)  is  obtained  as  the  difference  between  the  total 
nonlinear  term  (8)  and  the  resolved  nonlinear  term  (6). 

The  above  described  procedure  has  its  exact  counterpart  in  the  physical  space  represen¬ 
tation.  Inverse  Fourier  transform,  signified  by  tilde,  of  iVn(k)  (Eq.  (3))  is  the  sum  of  the 
convective  and  pressure  terms  in  the  Navier-Stokes  equation  in  the  physical  space  coordinates 


Nn(x)  =  -«j(x) 


diin(x)  dp(x) 


(9) 


dxi  dxn 

Similarly,  using  Nnuv( k)  we  can  define  its  physical  space  counterpart  N%v(x)  as  well  as 
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/Vf'V(x) 

which  is  the  inverse  Fourier  transform  of  .Vn^vVk)  truncated  to  the  region  K  (which 
is  either  C  or  S).  A  ^v(x)  can  be  interpreted  as  the  contribution  to  the  rate  of  change  of  the 
velocity  field  un(x)  at  a  point  x  made  by  the  nonlinear  interactions  involving  modes  from 
the  spectral  regions  U  and  V.  Note  that  these  interactions  influence  all  modes  k  which  can 
form  a  triangle  with  two  other  modes  such  that  one  is  in  U  and  the  other  in  V.  N^uv(x) 
represents  a  contribution  to  the  rate  of  change  of  un(x)  which  is  made  by  all  modes  from  K 
interacting  nonlinearly  with  modes  in  U  and  V.  Finally,  the  inverse  Fourier  transform  of  (5) 
is 


I«f(x)  =  ■W’fiffx)  +  ,v„(x|*c)  +  A'-fxIfc),  (10) 

where  the  resolved  nonlinear  term  IV„(x|fcc)  and  the  subgrid-scale  term  -/V*(x|fcc)  in  the 
physical  space  are  obtained  Fourier  transforming  (6)  and  (7),  respectively. 

In  the  LES  the  most  fundamental  requirement  is  that  the  models  employed  properly  ap¬ 
proximate  effects  of  subgrid-scale  interactions  on  the  energetics  of  the  resolved  scales.  Thus, 
in  assessing  the  models,  energy  equations  rather  than  momentum  equations  are  usually  con¬ 
sidered.  In  the  spectral  space  the  equation  for  the  energy  amplitude  ||u(k)|2  =  |un(k)u*(k) 
of  mode  k  obtained  from  (2)  is: 

|jMk)|2  =  -2^2t|u(k)|!  +  r(k),  (in 

where  T{ k)  is  the  nonlinear  energy  transfer 

T(k)  =  Re[un(k)Nn(k)}.  (12) 
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For  homogeneous  turbulence  the  above  equations  are  usually  considered  after  summing 
up  contributions  from  all  modes  with  a  prescribed  wavelength  | k j  =  k,  giving: 

-E(k)  =  -2vk2E(k)  +  T(k),  (13) 

where  E(k)  and  T(k)  are  the  classical  energy  and  transfer  spectra,  respectively,  for  homoge¬ 
neous  turbulence. 

Similarly,  the  detailed  energy  transfer  to/from  mode  k  caused  by  its  interactions  with 
wave  numbers  p  in  a  prescribed  region  V  of  the  wave  number  space  and  q  =  k  — p  in  another 
region  Q  is 

TVQ(k)  =  Re{un(k)NnVQ(k)}.  (14) 

Truncating  TPC(k)  to  another  region  K  results  in  the  quantity  TICVQ(k)  which  is  inter¬ 
preted  as  the  energy  transfer  to  the  region  1C  resulting  from  nonlinear  interactions  of  scales 
in  K  with  scales  in  V  and  Q.  For  homogeneous  turbulence  the  regions  V  and  Q  are  usually 
chosen  as  spherical  wave  number  bands  centered  at  wave  numbers  p  and  q,  respectively.  In 
this  case  quantity  TPC(k)  will  be  denoted  by  T(k\p,q).  Summing  quantity  T(k|p,  q)  over 
spherica'  shells  with  thickness  A k  =  1  centered  at  wave  number  k  gives  a  function  denotid 
by  either  TICVQ(k)  or  T(k\p,q) 

rCTe(i)  =  r(t|P,,)=  •£  r(k|P,,).  m 

k-%Ak<\k\<k+%&k 

Total  nonlinear  energy  transfer  T(k)  to  the  wave  number  band  k  is  obtained  by  summit;* 


contributions  T(k\p,  q)  from  ail  possible  bands  p  and  q : 

rm  =  ££!•(  %,,)  =  £/>(  Up).  (16) 

P  <3  P 

Here,  the  function  P{k\p)  is  a  result  of  summation  of  T(k\p,q)  over  all  g-bands  and  is 
interpreted  as  the  energy  transfer  between  wave  number  bands  k  and  p. 

With  this  notation,  the  energy  equation  for  the  energy  spectrum  Ec(k )  of  resolved  scales, 
obtained  from  Eq.  (5),  is 

jtEc{k)  =  -2  uk'1Ec{k)  +  T(k\kc)  +  T’{k\kc),  (17) 

where 

T(k\kc)  =  Tccc{k),  (18) 

and 

Ts{k\kc)  =  Tccs(k)  +  Tcss(k).  (19) 

Equivalent  expressions  in  the  physical  space  are  obtained  by  considering  an  equation  for 
the  rate  of  change  of  the  turbulent  energy  of  the  resolved  scales  e£(x)  =  |u£(x)tx£(x): 

=  ^(x)V^(x)  +  f  (xlfe)  +  7"(x|U  (20) 


where 


(21) 


T(x\kc)  =  u%{x)Nn(x\kc), 


is  the  resolved  energy  transfer  and 


f'Mfc)  =  if(x)A'*(x|te)  (22) 

is  the  subgrid-scale  energy  transfer  in  the  physical  space  representation. 

4  Results 

Using  the  methodology  described  in  the  previous  section  and  employing  the  numerically 
simulated  velocity  fields  it  is  possible  to  compute  directly  the  subgrid-scale  energy  transfer 
for  any  prescribed  cutoff  wave  number  kc  <  km.  It  is  customary  to  represent  spectral 
subgrid-scale  energy  transfer  in  terms  of  the  subgrid-scale  eddy  viscosity 


ue(k\kc) 


T’(k\kc ) 
2k2Ec(k)' 


k  <  kc, 


(23) 


which,  following  Kraichnan  is  usually  normalized  by  the  factor  equal  to  the  product  of 
the  velocity  scale  [f^fc^fc,.]1/2  and  the  length  scale  l/kc  at  the  cutoff  kc 


v  K(k\k  )  - 

e  Wkc)~  [E(kc)/kc]W 


(24) 


In  order  to  compute  the  function  T*(A:jfcc),  which  depends  on  the  length  of  the  wave 
number  k ,  according  to  Eq.  (15)  summation  over  all  wavevectors  k  in  a  thin  spherical  shell 
centered  at  k  must  be  performed.  The  components  of  such  a  sum  are  in  general  of  both 
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signs  implying  that  a  particular  mode  k  may  be  either  losing  energy  (forward  transfer)  or 
gaining  energy  (inverse  transfer)  because  of  the  subgrid-scale  nonlinear  interactions.  To 
assess  the  relative  importance  of  these  two  processes  we  have  performed  partial  summations 
over  components  of  same  sign,  effectively  splitting  the  subgrid-scale  energy  transfer  to/from 
scales  k  into  the  forward  and  the  inverse  transfer  contributions.  This  procedure  is  equivalent 
to  decomposing  the  total  eddy  viscosity  into  two  parts,  a  positive  one  associated  with  the 
forward  energy  transfer,  and  a  negative  one  associated  with  the  inverse  energy  transfer.  In 
Fig.  3  we  plot  the  total  eddy  viscosity  and  its  positive  and  negative  components  computed 
for  three  different  cutoff  wave  numbers:  one,  kc  =  20,  at  the  end  of  the  (nominal)  inertial 
subrange;  the  next,  kc  =  40,  at  the  beginning  of  the  dissipation  range;  and  the  last  one, 
kc  =  80,  deep  in  the  dissipation  range.  For  kc  =  20  the  total  eddy  viscosity  is  predominantly 
positive,  with  the  absolute  values  of  the  negative  component  about  30-50%  of  the  values 
of  the  positive  component  for  k/kc  <  0.6.  For  k/kc  >  0.6  the  ratio  of  the  negative  to  the 
positive  component  decreases  to  about  20%.  For  kc  =  40  in  the  range  k/kc  <  0.6  the  positive 
and  negative  components  nearly  balance  each  other  with  the  resulting  total  eddy  viscosity 
close  to  zero.  For  k/kc  >  0.6  both  components  exhibit  cusp-like  behavior  with  the  cusp  for 
the  positive  component  much  stronger  than  for  the  negative  one.  Nevertheless,  even  close 
to  the  cutoff  the  ratio  of  the  negative  to  the  positive  component  is  about  15%.  For  the  case 
kc  =  80  in  the  range  k/kc  <  0.6  the  positive  component  is  practically  zero  and  the  negative 
one  is  slightly  less  than  zero,  resulting  in  small  negative  values  of  the  total  eddy  viscosity. 
Beyond  that  range,  for  k  approaching  the  cutoff,  the  positive  component  increases  very 
rapidly,  reaching  at  the  cutoff  kc  values  by  factor  20  greater  than  the  values  of  the  negative 
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component.  In  the  last  figure  we  also  plot  the  subgrid-scale  eddy  viscosity  calculated  by 
Kraichnan  ^  and  Cholet  and  Lesieur  ^  from  the  analytical  theories  of  turbulence  under  the 
assumption  of  the  infinite  inertial  range.  This  function  is  essentially  constant  (equal  to  0.267) 
for  k/kc  <  0.6,  and  exhibits  the  cusp-like  behavior  for  k/kc  >  0.6.  We  conclude  from  this 
analysis  that  the  spectral  inverse  energy  transfer  may  be  quite  significant,  in  some  cases 
comparable  to  the  forward  transfer  for  given  scales  k.  However,  in  all  cases  the  forward 
transfer  dominates  as  the  cutoff  wave  number  is  approached.  Since  the  transfer  is  obtained 
by  multiplying  the  eddy  viscosity  by  k 2,  the  cusp  in  the  eddy  viscosity  for  k/kc  >  0.6  is 
actually  even  more  significant  for  the  subgrid-scale  transfer. 

Using  Eq.  (22)  we  have  computed  the  subgrid-scale  energy  transfer  in  the  physical  space 
T*(x|fcc)  for  several  spectral  cutoff  wave  numbers.  In  Fig.  4  we  plot  a  cross-section  of  this 
quantity  for  kc  =  20  and  kc  =  40  for  a  plane  in  the  impermeable  box  located  at  y  =  7r/4 
The  larger  spectral  cutoff  wave  number  results  in  presence  of  smaller  scales  in  the  physical 
space.  Regions  of  the  forward  transfer  (broken  contours)  and  the  inverse  transfer  (solid 
contours)  are  clearly  visible.  Even  though  the  overall  subgrid-scale  transfer  integrated  over 
the  computational  box  is  negative,  the  forward  and  inverse  transfer  regions  in  these  plots 
are  roughly  in  balance.  This  indicates  that  both  effects  may  be  equally  important  in  the 
dynamics  of  the  flow.  This  conclusion  agrees  with  the  corresponding  conclusion  reached  in 
the  analysis  of  the  spectral  subgrid-scale  transfer.  It  should  be  noted,  however,  that  there  is 
no  direct  relation  between  sets  of  spectral  modes  characterized  by  positive/negative  transfer 
and  the  physical  space  regions  with  the  same  characteristics.  Assuming  that  in  the  large 
scale  momentum  equation  (5)  the  subgrid-scale  nonlinear  term  Wn*(k|fcc)  is  represented  using 
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the  classical  spectral  eddy  viscosity  model  ue(k\kc)  of  Kraichnan  ^  and  Cholet  and  Lesieur  ^ 

Nna(k\kc)  =  - ve(k\kc)Pui(k ),  (25) 

we  have  calculated  the  physical  space  subgrid-scale  energy  transfer  for  this  model  from  (22). 
The  results  of  the  calculations  are  plotted  in  Fig.  5  for  kc  =  20  and  the  same  cross-plane  as 
in  Fig.  4.  The  modeled  transfer  is  predominantly  of  the  forward  type  as  expected  from  the 
use  of  the  strictly  positive  eddy  viscosity  but  the  appearance  of  weak  inverse  transfer  regions 
may  seem  surprising.  However,  it  should  be  noted  that  the  molecular  viscosity  term  in  the 
incompressible  Navier-Stokes  equations  results  in  two  distinct  effects  in  the  energy  equation: 
the  kinetic  energy  dissipation,  which  is  negative  everywhere,  and  the  change  in  the  kinetic 
energy  caused  by  work  done  by  viscous  stresses,  which  locally  in  space  may  be  either  positive 
or  negative.  Therefore  any  model  which  approximates  the  subgrid-scale  nonlinear  term  A’nJ 
by  a  viscous-like  term  in  the  Navier-Stokes  equations  may  contain  regions  of  the  increasing 
kinetic  energy  caused  by  work  done  by  the  modeled  stresses.  In  practice,  however,  as  seen  in 
Fig.  5,  these  positive  regions  are  quite  insignificant  since  they  occupy  much  less  space  than 
the  negative  regions  and  have  also  much  lower  maximum  values.  Obviously,  such  models 
give  poor  representation  of  the  actual  subgrid-scale  energy  transfer  as  seen  comparing  tin- 
actual  and  modeled  transfers  shown  in  Fig.  4(a)  and  Fig.  5,  respectively.  The  conclusion-* 
from  the  physical  space  analysis  of  the  subgrid-scale  energy  transfer  parallel  those  drawn 
from  the  spectral  space  analysis:  a  relative  importance  of  the  inverse  energy  transfer  procr-*- 
and  an  inability  of  the  classical  subgrid-scale  models  to  properly  account  for  it. 

Cusps  observed  in  spectral  eddy  viscosities  in  the  vicinity  of  the  cutoff  wave  numb.-r 
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suggest  that  the  total  energy  transfer  across  this  wave  number  is  dominated  by  energy 
exchanges  among  resolved  and  unresolved  scales  from  the  vicinity  of  the  cutoff.  Indeed, 
it  has  been  established  in  a  number  of  papers  that  in  numerically  simulated 

turbulence  at  low  Reynolds  numbers  the  energy  transfer  beyond  the  energy  containing  range 
is  local,  occuring  between  scales  of  similar  size,  even  though  the  nonlocal  wave  number 
triads  with  one  scale  in  the  energy  containing  range  are  responsible  for  this  local  transfer. 
One  would  thus  expect  that  the  subgrid-scale  nonlinear  interactions  between  the  resolved 
scales  ( k  <  kc)  and  the  unresolved  scales  characterized  by  wave  numbers  slightly  greater 
than  the  cutoff  wave  number  kc  will  dominate  the  subgrid-scale  energy  transfer  process.  To 
evaluate  this  hypothesis  in  more  detail  we  have  calculated,  for  several  values  of  the  cutoff 
wave  number  kc,  the  subgrid-scale  energy  transfer  for  the  truncated  velocity  fields  obtained 
from  the  original  field  by  setting  to  zero  all  modes  with  wave  numbers  k  >  ckc,  where  c  was 
equal  to  3/2  and  2.  In  this  way  the  effect  of  all  modes  k  >  ckc  on  the  subgrid-scale  energy 
transfer  is  eliminated.  In  Fig.  6  we  plot  the  resulting  spectral  subgrid-scale  eddy  viscosities 
(24)  for  kc  =  20  and  kc  =  40  and  compare  them  with  the  eddy  viscosities  computed  using 
the  full  velocity  field,  i.e.  with  all  modes  k  <  km  being  non-zero.  It  is  seen  that  the  value 
of  the  eddy  viscosity  computed  for  c  =  3/2  provides  a  very  good  approximation  to  the  total 
eddy  viscosity  while  for  c  =  2  both  quantities  are  practically  indistinguishable  on  the  plots. 
The  similarly  calculated  subgrid-scale  transfer  in  the  physical  space  (22)  is  shown  in  Fig. 
7  for  the  cutoff  wave  number  kc  =  20  and  two  values  of  the  parameter  c,  3/2  and  2.  The 
plane  shown  is  the  same  as  in  Fig.  4.  The  spatial  structure  of  the  subgrid-scale  energy 
transfer  in  Fig.  7  and  Fig.  4(a)  is  the  same,  with  differences  seen  only  in  the  values  of  the 


18 


transfer  at  particular  locations.  For  c  =  3/2  the  peak  values  of  the  approximated  transfer 
(Fig.  7(a))  may  depart  by  about  10%  from  the  exact  values  (Fig.  4(a))  with  the  departures 
decreasing  to  about  5%  for  c  =  2  (Fig.  7(b)).  Therefore,  both  in  the  spectral  and  the 
physical  space  representation  the  subgrid-scale  energy  transfer  for  the  resolved  modes  k  <  kc 
can  be  determined  with  high  accuracy  by  considering  their  interactions  with  a  limited  range 
of  unresolved  modes  kc  <  k  <  2 kc.  It  may  be  instructive  to  note  that  that  for  kc  =  20  and 
km  =  170  the  resolved  modes  constitute  about  0.0016  of  all  modes,  and  modes  with  k  <  2 kc 
about  0.013  of  all  modes.  Thus  the  dynamics  of  the  largest  0.16%  modes  is  determined 
almost  entirely  by  their  nonlinear  interactions  with  about  1%  of  all  modes,  the  remaining 
99%  modes  not  affecting  visibly  the  largest  scales.  Moreover,  the  lack  of  direct  influence  of 
small  scales  k  >  2 kc  on  the  energetics  of  the  large  resolved  scales  k  <  kc  implies  that  the 
direct  nonlocal  energy  transfer,  inherent  in  the  classical  eddy  viscosity  theories  ,  js  not 
present  in  our  simulations.  The  dynamics  of  the  largest  modes  observed  in  the  simulations  is 
quite  similar  to  the  classical  picture  of  the  dynamics  of  the  energy  containing  range  in  high 
Reynolds  number  turbulence.  Quoting  Batchelor  ^ :  ”It  seems  that  the  energy-containing 
eddies  determine  the  rate  of  energy  transfer  by  their  mutual  interactions,  and  the  larger 
wave-numbers  adjust  themselves,  according  to  the  Reynolds  number,  in  order  to  convert 
this  energy  into  heat  at  the  required  rate.” 

For  the  purpose  of  subgrid-scale  modeling  it  is  important  to  investigate  relations  between 
observed  subgrid-scale  energy  transfer  and  various  features  of  the  resolved  velocity  field 
(4).  In  Figs.  8,  9,  and  10  we  plot  in  the  physical  space  representation  the  kinetic  energy 
|u(x)  •  u(x),  the  enstrophy  \u>(x)  •  u>(x),  where  u;(x)  is  the  vorticity,  and  the  dissipation 


19 


rate  ^(du./chr*  +  duk/dij)2,  respectively.  All  these  quantities  are  computed  using  both 
the  full  velocity  field  (all  modes  k  <  km  are  nonzero)  and  the  velocity  field  (4)  truncated 
at  kc  =  20.  One  cross-sectional  plane  at  the  location  y  =  (3/4)7r  in  the  impermeable  box 
is  plotted  and  for  comparison  we  also  plot  in  Fig.  11  the  subgrid-scale  energy  transfer  in 
the  same  plane.  Spatial  structure  of  energy  fields  for  the  full  and  truncated  fields  is  nearly 
the  same.  This  feature  is  expected  since  the  modes  k  <  kc  contain  most  of  the  total  energy. 
The  peak  values  of  the  energy  for  the  full  field  may  exceed  by  30%  the  peak  values  for  the 
truncated  field.  The  maximum  values  for  the  dissipation  and  the  enstrophy  fields  computed 
using  the  full  velocity  fields  are  by  a  factor  5  greater  than  for  the  truncated  fields,  indicating 
fairly  large  contributions  coming  from  higher  wave  numbers  k  >  kc.  The  importance  of 
these  wave  numbers  is  also  reflected  in  the  spatial  structure  of  these  quantities,  with  the  full 
fields  showing  the  presence  of  much  smaller  scales  than  the  truncated  fields.  Despite  these 
differences  between  the  full  and  the  truncated  fields,  for  the  enstrophy  the  spatial  structure 
of  the  large  scale  component  (Fig.  9(b))  is  remarkably  similar  to  the  structure  of  the  total 
enstrophy  (Fig.  9(a)).  In  particular  the  regions  of  large  values  of  the  total  enstrophy  are 
very  well  correlated  with  the  regions  where  the  truncated  field  also  gives  large  values.  This 
result  is  somewhat  surprising  since  the  large  scale  enstrophy  field  is  determined  using  only 
0.16%  of  all  modes.  It  suggests  that  these  largest  scales  contain  most  of  phase  information 
required  to  determine  spatial  structure  of  the  enstrophy  field,  and  the  role  of  higher  wave 
number  modes  is  to  merely  reflect  the  fact  that  the  velocity  gradients  are  steeper  than  can 
be  resolved  by  the  low  wave  number  modes.  In  other  words,  larger  wave  numbers  in  the 
spectral  space  are  needed  to  resolve  steep  velocity  gradients  rather  them  small  eddies  thought 
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of  as  small,  individual  flow  structures  like  localized  vortices.  The  level  of  correlation  between 
the  full  and  truncated  fields  for  the  dissipation  is  lower.  The  regions  of  the  most  intense 
dissipation  for  the  truncated  field  (located  along  diagonals,  half  way  between  the  center  and 
the  corners  of  the  plotted  plane  in  Fig.  10(b))  correlate  well  with  the  full  dissipation  field 
in  the  same  region  (Fig.  10(a))  but  some  equally  strong  regions  in  the  full  field  farther 
away  from  the  center  do  not  have  clear  counterparts  in  the  truncated  field.  Finally,  all  three 
quantitites  computed  using  the  truncated  velocity  field  were  compared  with  the  subgrid- 
scale  energy  transfer  plotted  in  Fig.  11.  There  is  some  level  of  spatial  correlation  between 
the  subgrid-scale  transfer  and  the  enstrophy  and  the  dissipation  fields,  with  the  regions  of 
significant  transfer  in  the  vicinity  (but  not  on  the  top  of)  regions  of  large  enstrophy  and 
dissipation.  Also,  the  regions  of  intense  large  scale  dissipation  are  usually  located  on  ihe 
peripheries  of  the  regions  of  intense  large  scale  enstrophy.  Interestingly,  the  regions  of  large 
subgrid-scale  transfer  seem  to  correlate  best  with  the  regions  of  large  scale  energy  (Fig.  9(b)). 
Such  correlations  were  observed  previously  ^  for  different  velocity  fields  but  no  convincing 
physical  explanation  of  this  observation  is  known.  The  above  observations  are  based  on  visual 
inspection  of  countour  plots  and  thus  have  a  very  qualitative  character.  A  more  quantitative 
procedure  would  have  to  be  used  to  evaluate  correlations  in  a  systematic  way.  Nevertheless, 
this  qualitative  analysis  clearly  illustrates  a  fairly  complex  character  of  inter-relations  among 
different  physical  quantities  and  gives  no  indications  that  any  simple  expression  for  the 
subgrid-scale  transfer  in  terms  of  the  resolved  energy,  enstrophy,  or  dissipation  exists.  We 
conclude  from  the  analysis  of  the  truncated  fields  that  the  the  subgrid-scale  energy  transfer 
is  at  best  marginally  correlated  with  the  large  scale  energy,  enstrophy,  and  dissipation.  This 
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analysis  also  reveals  that  the  turbulent  activity  is  spatially  intermittent  and  its  physical 
locations  are  determined  by  the  mutual  interactions  of  the  largest  scales.  The  high  wave 
number  modes  in  this  flow  cannot  be  interpreted  as  individual,  small  scale  turbulent  eddies, 
but  reflect  the  presence  of  steep  gradients  at  the  spatial  locations  determined  by  the  large 
scales. 

The  observed  importance  of  the  large  scales,  which  constitute  only  a  minute  fraction 
of  all  spectral  modes,  in  the  dynamics  of  turbulence  is  encouraging  since  it  suggests  that 
their  dynamics  may  be  almost  self-contained  and  thus  the  accurate  subgrid-scale  models 
based  on  the  large  scale  velocity  information  should  be  possible.  The  term  "almost  self- 
contained  dynamics”  is  not  very  precise  but  can  be  illustrated  by  the  following  example.  In 
Fig.  12(a)  we  plot  one  plane  from  the  resolved  nonlinear  transfer  field  (Eq.  (21))  and  in 
Fig.  12(b)  the  corresponding  result  for  the  total  nonlinear  transfer,  i.e.  the  sum  of  (21) 
and  (22).  It  is  seen  that  the  spatial  structure  of  the  resolved  nonlinear  transfer  is  highly 
correlated  with  the  structure  of  the  total  nonlinear  transfer.  In  that  sense  the  dynamics  of 
the  large  resolved  scales,  which  involves  interactions  with  all  modes  in  the  system  (the  total 
transfer)  is  "almost”  the  same  as  the  internal  nonlinear  dynamics  of  the  large  scales  on  I  > 
(the  resolved  transfer).  The  difference  between  both  quantities  is  of  course  the  subgrid-scale 
nonlinear  transfer  (22),  which  when  viewed  this  way,  is  a  small  correction  to  the  resolved 
transfer  needed  to  get  the  total  transfer  and  to  account  for  the  non-conservative  character 
of  the  entire  system.  The  large  degree  of  correlation  between  the  resolved  and  the  tot  a; 
transfers  suggests  that  in  the  subgrid-scale  modeling  a  good  strategy  may  be  to  model  the 
total  (unknown)  transfer  using  the  resolved  (known)  transfer  rather  than  to  model  their 


difference,  the  subgrid-scale  transfer,  as  is  always  done. 


5  Conclusions 

We  have  performed  a  detailed  analysis  of  the  effects  of  the  subgrid-scale  nonlinear  interactions 
on  the  energetics  of  isotropic  turbulence.  The  analyzed  turbulent  velocity  field  was  obtained 
from  a  direct  numerical  simulation  of  the  Taylor-Green  vortex  flow.  Symmetries  of  the  flow 
allowed  to  reach  the  spatial  resolution  in  the  simulation  equivalent  to  5123  mesh  points  and 
the  Reynolds  number  R\  ~  70.  At  this  Reynolds  number  the  flow  exhibits  a  beginning  of 
the  inertial  range  dynamics  at  the  lowest  wave  numbers.  However,  even  these  low  as  well  as 
all  higher  wave  numbers  are  still  dominated  by  dissipative  processes.  Therefore,  while  our 
conclusions  are  certainly  valid  for  the  disspation  range  dynamics  it  is  less  certain  that  they 
are  applicable  to  the  inertial  range  dynamics. 

An  important  feature  of  the  computed  subgrid-scale  energy  transfer,  in  both  spectral 
and  physical  space  representation,  is  the  presence  of  significant  inverse  energy  transfers, 
from  the  unresolved  to  the  resolved  scales.  The  inverse  subgrid-scale  transfer  was  predicted 
and  observed  before  in  the  context  of  spectral  dynamics  of  homogeneous  turbulence  4»6,11,8 
and  in  the  physical  space  for  homogeneous  and  wall-bounded  turbulent  flows  The 

observed  significance  of  the  inverse  transfer  in  the  energetics  of  the  resolved  scales  implies 
that  successful  subgrid-scale  models  should  properly  account  for  such  effects.  At  the  present 
time  these  effects  are  rarely  taken  into  account  in  the  subgrid-scale  modeling  procedures. 
If  accounted  for  they  are  modeled  by  either  adding  a  random  force  to  the  subgrid-scale 
equations  or  extrapolating  from  the  dynamics  of  the  resolved  scales  Since  most  of  the 
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subgrid-scale  transfer  observed  in  this  work  is  caused  by  interactions  among  highly  correlated 
modes  on  both  sides  of  the  cutoff  wave  number,  approximating  effects  of  such  interactions  by 
random  forces  is  debatable.  An  approach  used  in  the  dynamic  subgrid-scale  model  ®  seems 
more  appropriate  but  it  suffers  from  modeling  the  inverse  transfer  by  a  diffusion  type  term 
with  a  negative  diffusion  coefficient,  mathematically  an  inherently  unstable  situation.  It 
appears  that  alternate  ways  of  the  subgrid-scale  modeling  which  overcome  these  conceptual 
and  mathematical  difficulties  should  be  explored. 

Our  analysis  also  reveals  that  the  nonlinar  dynamics  of  the  resolved  modes  with  wave 
numbers  k  <  kc  is  governed  almost  exclusively  by  their  interactions  with  a  limited  range 
of  modes  with  wave  numbers  not  exceeding  2 kc  and  nonlocal,  eddy-viscosity  type  energy 
transfer  is  not  observed.  Thus,  in  agreement  with  the  classical  picture  of  the  turbulence 
dynamics  the  large  scales  o'"  a  turbulent  flow  determine  the  energy  flux  down  the  spectrum 
and  the  small  scales  play  entirely  passive  role  by  adjusting  themselves  in  such  a  way  as  to 
accomodate  this  energy  flux  prescribed  by  the  large  scales. 

The  physical  space  energy,  enstrophy,  and  dissipation  have  been  computed  for  the  full  and 
truncated  velocity  fields  and  compared  with  the  subgrid-scale  energy  transfer  for  the  same 
truncation  wave  number.  Surprisingly,  these  physical  quantities  computed  for  both  full  and 
truncated  fields  show  many  similar  spatial  features  despite  the  fact  that  the  truncated  field 
contains  only  0.16%  of  all  modes  present  in  the  system.  This  result  reinforces  our  conclusion 
about  the  dominant  role  played  by  the  very  largest  scales  in  the  dynamics  of  the  flow.  The 
level  of  correlation  between  these  quantities  and  the  subgrid-scale  transfer  varies  from  weak 
for  the  enstrophy  and  dissipation,  to  moderate  for  the  energy. 
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Finally,  using  the  observed  importance  of  the  nonlinear  interactions  among  the  largest 
scales  in  the  overall  dynamics  of  the  resolved  scales  it  is  i  iggested  that  the  resolved  nonlinear 
term  may  possibly  serve  as  a  basis  of  a  new  approach  to  the  subgrid-scale  modeling. 
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Figure  Captions 

Figure  1.  Time  evolution  of  flow  quantities:  (a)  total  kinetic  energy;  (b)  total  dissipation; 
(c)  skewness;  (d)  microscale  Reynolds  number. 

Figure  2.  Spectral  quantities  at  the  end  ot  the  run  ( t  =  IS):  (a)  unnormalized  spectra  of 
the  energy  (solid  line)  and  the  dissipation  (broken  line);  (b)  the  dissipation  spectrum  (solid 
line)  and  the  transfer  spectrum  (broken  line)  outside  the  energy  containing  range.  Both 
quantities  plotted  using  a  log-linear  scale  to  accentuate  their  exponential  behavior. 

Figure  3.  The  spectral  subgrid-scale  eddy  viscosity  (solid  line)  and  its  negative  (broken 
line)  and  positive  (dotted  line)  components:  (a)  kc  =  20;  (b)  kc  =  40;  (c)  kc  =  80. 

Figure  4.  The  subgrid-scale  energy  transfer  in  the  physical  space  representation:  (a) 
kc  =  20;  (b)  kc  =  40.  Plane  y  =  7r/4  in  the  impermeable  box  is  shown.  Here  and  in 
all  subsequent  contour  plots  the  solid  lines  represent  positive  values  and  the  broken  lines 
represent  negative  values. 

Figure  5.  The  subgrid-scale  energy  transfer  in  the  physical  space  computed  using  the 

4  .  5 

spectral  eddy  viscosity  model  of  Kraichnan  and  Chollet  and  Lesieur  . 

Figure  6.  The  spectral  subgrid-scale  eddy  viscosity:  (a)  kc  =  20;  (b)  kc  =  40.  The 
velocity  fields  used  to  compute  this  quantity  were  the  full  field  (solid  line),  the  full  field 
truncated  at  (3/2 )kc  (broken  line),  and  the  full  field  truncated  at  2 kc  (dotted  line). 

Figure  7.  The  physical  space  subgrid-scale  energy  transfer  computed  for  kc  =  20  and  the 
full  fields  truncated  at  (a)  (3/2 )kc  and  (b)  2 kc. 

Figure  8.  The  kinetic  energy  field  in  a  plane  y  =  (3/4)7 r  in  the  impermeable  box:  (a) 
computed  using  the  full  velocity  field;  (b)  computed  using  the  resolved  velocity  field  with 
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the  truncation  wave  number  kc  =  20. 


Figure  9.  The  enstrophy  field  in  a  plane  y  =  (3/4)7r  in  the  impermeable  box:  (a) 
computed  using  the  full  velocity  field;  (b)  computed  using  the  resolved  velocity  field  with 
the  truncation  wave  number  kc  =  20. 

Figure  10.  The  dissipation  field  in  a  plane  y  =  (3/4)t  in  the  impermeable  box:  (a) 
computed  using  the  full  velocity  field;  (b)  computed  using  the  resolved  velocity  field  with 
the  truncation  wave  number  kc  =  20. 

Figure  11.  The  subgrid-scale  energy  transfer  in  the  physical  space  representation  com¬ 
puted  for  the  cutoff  wavenumber  kc  =  20  and  shown  in  a  plane  y  =  (3/4)7r  in  the  impermeable 
box. 

Figure  12.  The  nonlinear  energy  transfer  to/from  the  resolved  modes  k  <  20  represented 
in  the  physical  space:  (a)  caused  by  interactions  with  the  resolved  modes  only;  (b)  caused 
by  interactions  with  the  resolved  and  unresolved  modes. 
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Abstract 


Using  results  of  direct  numerical  simulations  of  isotropic  turbulence  the  subgrid-scale  energy 
transfer  in  the  physical  space  is  calculated  exactly  employing  a  spectral  decomposition  of  the 
velocity  field  into  large  and  small  scales.  Both  the  forward  and  the  inverse  subgrid-scale 
transfer  components  are  found  to  be  significant.  The  spatial  structure  of  the  exact  subgrid- 
scale  transfer  is  qualitatively  compared  with  the  spatial  structure  of  a  number  of  physical 
quantities  which  are  considered  to  govern  the  dynamics  of  the  large  scales  of  turbulence.  We 
find  that  all  quantities  determined  by  the  first  derivatives  of  the  velocity  field  correlate  poorly 
with  the  transfer  which  is  largest  at  the  peripheries  of  regions  characterized  by  large  values 
of  these  quantities.  The  spatial  structure  of  the  transfer  correlates  much  better  with  the 
large  scale  energy  and  the  Smagorinsky’s  subgrid-scale  energy  transfer  which  is  determined 
by  the  second  derivatives  of  the  velocity  field.  None  of  the  considered  quantities  is  capable  of 
predicting  sign  of  the  subgrid-scale  transfer. 


1  Introduction 


It  is  generally  recognized  that  in  the  forseeable  future  numerical  predictions  of  turbulent  flows  at 
high  Reynolds  numbers  will  have  to  rely  on  models  accounting  for  the  effects  of  nonlinear  interactions 
between  resolved  (large)  and  unresolved  (small)  scales  of  motion.  Commonly  used  models,  either 
for  Reynolds  averaged  Navier-Stokes  (RANS)  equations  or  for  large  eddy  simulations  (LES),  are 
known  to  be  deficient  in  many  respects  and  the  improvements  in  the  models  will  be  difficult  to 
achieve  in  the  absence  of  underlying,  generally  accepted  theory  of  turbulence.  This  gives  impetus 
to  fundamental  investigations  of  nonlinear  interactions  in  turbulent  flows  which  may  eventually 
provide  better  models.  In  such  investigations  the  nonlinear  interactions  between  prescribed 
scales  of  motion  are  analyzed  using  accurately  resolved  velocity  fields  obtained  in  direct  numerical 
simulations.  In  the  context  of  large  eddy  simulations  a  class  of  subgrid-scale  nonlinear  interactions 
is  of  prime  importance.  Clark  et  al.  ^  pioneered  use  of  exactly  computed  subgrid-scale  interactions 
to  assess  phenomenological  subgrid-scale  eddy  viscosity  models  in  the  physical  space  representation. 

t 

Later,  Domaradzki  et  al.  °  applied  similar  methodology  to  compare  the  exactly  computed  spectral 
subgrid-scale  eddy  viscosity  with  eddy  viscosities  predicted  by  the  analytical  theories  of  turbulence. 
More  recently  Domaradzki  et  al.  D  used  the  same  approach  in  investigating  the  inverse  subgrid-scale 
energy  transfer  in  both  spectral  and  physical  space  representation  and  the  degree  of  localness  of  the 
subgrid-scale  transfer. 

Apart  from  the  usefulness  of  this  approach  in  assessing  existing  subgrid-scale  models  we  be¬ 
lieve  that  it  may  also  be  very  useful  in  investigation  of  coherent  structures.  Indeed,  despite  years 
of  research  devoted  to  coherent  structures  in  turbulent  flows  controversies  persist  a s  to  their  dy- 
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namical  significance  in  ('volution  of  turbulence  with  the  resulting  lack  ot  progress  in  <i<". I'iiunnii 

turbulence  models  that  would  explicitly  include  information  about  coherent  structures.  Since  the 

nonlinear  energy  transfer  is  a  principal  physical  process  influencing  the  evolution  of  turbulent  fields 

the  dynamical  significance  of  the  coherent  structures  may  be  established  if  it  is  shown  that  they  are 

associated  with  significant  transfer.  More  generally,  one  may  argue  that  the  existence  of  a  noticable 

correlation  between  a  given  physical  quantity  and  the  transfer  implies  the  dynamical  significance  of 

£* 

this  quantity.  Along  these  lines  Domaradzki  et  al.  D  investigated  qualitative  correlations  between 
subgrid-scale  energy  transfer  and  the  energy,  vorticity,  and  dissipation  fields  for  the  symmetric 
Taylor-Green  vortex  flow. 

In  this  paper  the  same  approach  is  used  to  investigate  the  subgrid-scale  energy  transfer  in  regular, 
non-symmetric  flows,  and  its  spatial  correlations  with  several  physical  quantities  which  customarily 
are  considered  as  being  dynamically  important  in  turbulent  flows. 

2  Numerical  simulations 

Description  of  field  A  -  1283. 

Description  of  field  B  -  2563. 

3  Basic  quantities 

Details  of  calculations  of  the  subgrid-scale  interactions  from  the  results  of  direct  numerical  simula¬ 
tions  of  isotropic  turbulence  are  described  by  Domaradzki  et  al.  Here  we  quote  only  the  main 
formulas.  For  homogeneous  turbulence  incompressible  Navier-Stokes  equations  in  spectral  (Fourier) 
representation  are: 
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•^-un(k)  =  —  j/fc2tt„(k)  +  iVn(k). 


(1) 


Here,  w„(k)  is  the  velocity  field  in  spectral  space,  with  the  explicit  dependence  on  time  omitted,  v 
is  the  kinematic  viscosity,  and  Nn( k)  is  the  nonlinear  term 


A'n(k)  =  -^F„/m(k)  J  rf3pu,(p)um(k  -  p)  ,  (2) 

where  tensor  P„/m(k)  accounts  for  the  pressure  and  incompressibility  effects.  The  summation  con¬ 
vention  is  assumed  throughout. 

The  wave  number  space  is  divided  into  two  non-overlapping  regions,  £  (|k|  <  kc)  signifying 
large,  resolved  scales,  and  S  (|kj  >  kc)  signifying  small,  unresolved  scales.  An  evolution  equation 
for  the  resolved  scales  i.e.  the  velocity  field  u„(k)  truncated  to  the  region  £ 


is: 


«f(k) 


un(k)  if  ke£ 

0  otherwise 


(3) 


|u/(k)  =  -uk\nc{  k)  +  Nn(k\kc)  +  Nn’(k\kc),  (4) 

where  jVn(k|fcc)  is  the  resolved  nonlinear  term  and  A^„*(k|^c)  is  the  subgrid-scale  nonlinear  term. 
These  terms  are  computed  as  follows.  First,  Eq.  (2)  is  used  with  the  full  (untruncated;  velocity 
fields  u/(p)  and  um(k  —  p)  and  the  result  is  truncated  to  the  region  £  to  obtain  the  total  nonlinear 
term 
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iVn‘°‘(k|Jbc)  =  Nn(k\kc)  +  Nns(k\kc) 


Next,  Eq.  (2)  is  used  again  with  the  truncated  velocity  fields  uf  (p)  and  u^(k  —  p)  and  the  result  is 
truncated  to  the  region  C  giving  the  resolved  nonlinear  term  iVn(k|A:c).  The  subgrid-scale  nonlinear 
term  is  obtained  as  the  difference  between  the  total  nonlinear  term  (5)  and  the  resolved  nonlinear 
term. 

The  above  described  procedure  has  its  exact  counterpart  in  the  physical  space  representation. 
Inverse  Fourier  transform,  signified  by  tilde,  of  JV„(k)  (Eq.  (2))  is  the  sum  of  the  convective  and 
pressure  terms  in  the  Navier-Stokes  equation  in  the  physical  space  coordinates 


iVn(x)  =  -u,-(x) 


dun{x)  dp(x) 
dxi  dxn 


Similarly,  the  inverse  Fourier  transform  of  (4)  is 


| =  i/V2uf  (x)  +  Nn(x\kc)  +  N^(x\kc),  (7) 

where  the  resolved  nonlinear  term  Nn(x\kc)  and  the  subgrid-scale  term  N*{x\kc)  in  the  physical 
space  are  obtained  Fourier  transforming  corresponding  expressions  in  the  spectral  space.  An  equa¬ 
tion  for  the  rate  of  change  of  the  turbulent  energy  of  the  resolved  scales 


E(x)  =  ^(x)5f(x) 

is  obtained  from  (7)  as 


(8) 
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(9) 


dE(x) 

dt 


vucn(x)V2rt(x)  +  f(x\kc)  +  fs(x\kc), 


where 


T{x\kc)  =  u*(x)Nn{x\ke), 


is  the  resolved  energy  transfer  and 


(10) 


r'MW  =  sf(x)jv;(x|ie)  (ii) 

is  the  subgrid-scale  energy  transfer  in  the  physical  space  representation. 

The  goal  of  the  subgrid-scale  modeling  is  to  obtain  as  good  as  possible  an  approximation  to  the 
subgrid-scale  nonlinear  term  7V*(x|fcc)  using  only  information  available  in  the  resolved  scales.  Thus 
on  the  level  of  the  energy  equation  we  are  interested  in  relations  between  observed  subgrid-scale 
energy  transfer  (11)  and  various  physical  quantities  computed  for  the  resolved  velocity  field  (3).  In 
addition  to  the  large  scale  energy  (8)  the  following  large  scale  quantities  have  also  been  considered: 
enstrophy,  dissipation,  vorticity  production,  enstrophy  production,  pressure,  the  second  invariant, 
and  Smagorinsky’s  subgrid-scale  energy  transfer.  These  quantities  are  defined  as  follows.  The 
enstrophy  is 


O(x)  =  -ft,(x)ft,(x)  , 

where  ft,  is  the  large  scale  vorticity.  Introducing  the  large  scale  rate-of-strain  tensor 


(12) 
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(13) 


Sij  =  duf/dij  +  duf/dxt) 


and  the  large  scale  rotation  tensor 


Rij  =  \{dufldxj  -  duj/dxi), 

we  define  dissipation 


(14) 


D{x)  =  ^5,j(x)5jj(x),  (15) 

the  square  of  the  vorticity  production  (i.e.  the  square  of  the  vortex  stretching  term  in  the  vorticity 
equation) 


P(x)  =  Q,(x)S,Jx)Slt(x)ih(x) , 


(16) 


enstrophy  production 


V(x)  =  a(x)5,J(x)0J(x) , 


(17) 


the  second  invariant  of  the  large  scale  velocity  gradient  tensor 


R(x)=l-lRi,R,j-S,,S„),  (18) 

and  the  quantity  proportional  to  the  subgrid-scale  energy  transfer  predicted  by  the  Smagorinsky’s 
model  ^ 
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4  Results 

Structures  in  isotropic  turbulence  are  often  identified  as  regions  of  concentrated  vorticity.  Visual- 
izations  of  numerically  simulated  isotropic  turbulence  ' show  existence  of  vortex  tubes  and 
vortex  sheets  in  such  flows  with  vortex  tubes  being  predominant.  In  Fig.  1  we  plot  isocontours 
of  enstrophy  for  the  field  A  [vorticity  level  is  ?].  The  vorticity  indeed  forms  elongated,  tube-like 
structures,  in  accordance  with  other  simulations  of  isotropic  turbulence.  In  Figs.  2(a)  and  2(b) 
we  plot  cross-sections  through  this  enstrophy  field  taken  in  horizontal  planes  at  distances  from  the 
lower  plane  equal  to  a  quarter  and  a  half  of  the  periodicity  length,  respectively.  In  Fig.  2(a)  two 
regions  of  intense  vorticity  are  clearly  visible.  The  elongated  region  in  the  lower  left  corner  results 
from  the  cut  of  the  vortex  tube  along  its  horizontally  oriented  axis  and  the  more  circular  region  to 
the  right  results  from  the  cut  of  the  vortex  tube  with  its  axis  inclined  with  respect  to  the  horizon¬ 
tal  plane.  In  Fig.  2(b)  several  regions  with  large  vorticity  values  have  generally  oval  shapes  and, 
after  comparison  with  Fig.  1  is  made,  appear  to  result  from  cutting  vortex  tubes  with  axes  in  the 
directions  not  far  from  normal  to  the  cutting  plane. 

In  Figs.  3(a)  and  3(b)  we  plot  subgrid-scale  energy  transfer  (11)  in  horizontal  planes  at  the  same 
locations  as  used  in  plotting  the  enstrophy.  An  important  observation  is  that  the  calculated  subgrid- 
scale  transfer  contains  regions  of  large  forward  transfer  (broken  line  contours)  and  regions  of  inverse 
transfer  of  comparable  magnitude  (solid  line  contours).  This  observation  is  consistent  with  results 
of  Piomelli  et  al.  11  for  turbulent  channel  flow  and  results  of  Domaradzki  et  al.  D  for  the  turbulent 
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Tay lor- Green  vortex  How.  The  presence  of  significant  inverse  transfer  is  in  direct  contradiction  'u 
the  implicit  assumptions  made  in  turbulence  models  that  the  subgrid-scale  interactions  are  of  purely 
dissipative  character  i.e.  transferring  energy  always  in  one  direction,  from  the  large  to  the  small 
scales. 

In  this  work  we  are  mainly  interested  in  how  regions  of  the  large  transfer  are  located  with  respect 
to  regions  characterized  by  by  large  values  of  the  other  physical  quantities.  In  order  to  investigate 
such  spatial  correlations  between  the  computed  subgrid-scale  energy  transfer  and  the  dynamics 
of  large  scales  we  plot  transfer  superimposed  on  contour  plots  of  the  previously  defined  functions 
of  the  large  scale  velocity  field.  No  attempt  was  made  to  make  point-wise  comparison  between 
the  exact  values  for  the  plotted  quantities  nor  the  point-wise  correlation  coefficient  was  computed. 
Our  analysis  is  qualitative  but  in  view  of  the  complexity  of  the  analyzed  fields  is  necessary  to 
establish  guidelines  for  more  quantitative  analysis  which  should  be  attempted  later.  For  instance, 
the  presence  of  both  positive  and  negative  values  for  the  computed  transfer  may  imply  negligible 
correlation  coefficient  between  the  transfer  and  a  purely  positive  quantity,  even  though  both  may 
have  a  very  similar  spatial  structure.  The  correlation  coefficient  will  also  be  small  if  one  quantity 
is  concentrated  in  the  vicinity  of  the  other  but  not  on  top  of  it.  Yet  such  relation  may  imply 
that  both  quantities  are  dynamically  correlated  even  though  the  pointwise  correlation  coefficient  is 
small.  Therefore  it  appears  that  the  qualitative  analysis  reported  here  is  a  reasonable  first  step  in 
establishing  approximate  correlations  between  various  fields  and  in  suggesting  directions  for  a  more 
quantitative  analysis. 

In  Fig.  4  large  values  of  enstrophy  (12)  are  usually  accompanied  by  large  positive/negative 
values  of  the  subgrid-scale  transfer  but  there  is  only  partial  overlap  between  these  quantities.  It 
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appears  that  maxima  of  both  fields  are  shifted  with  respect  to  '*ach  other  Midi  t  iiat  ivsiioU'  .<f 
transfer  are  located  on  the  edges  of  high  enstrophy  regions.  Even  though  transfer  seems  to  be 
associated  with  enstrophy  the  reverse  is  not  necessarily  true  i.e.  the  regions  of  high  transfer  may 
be  present  without  any  significant  enstrophy  in  their  vicinity. 

Higher  degree  of  correlation  is  observed  in  Fig.  5  between  the  subgrid-scale  transfer  and  the  large 
scale  energy  (8).  In  this  case  the  regions  of  the  most  intense  subgrid-scale  transfer  lie  almost  directly 
on  top  of  high  energy  regions.  This  is  particularly  clearly  illustrated  in  Fig.  5(b)  where  a  sequence 
of  the  subgrid-scale  regions  of  alternating  signs  forms  a  pattern  overlapping  with  the  pattern  formed 
by  the  large  scale  energy.  The  observed  correlation  is  puzzling  since  it  is  expected  that  the  energy 
transfer  process  is  associated  with  gradients  of  the  velocity  field,  not  with  the  velocity  field  itself. 
This  correlation  may  be  explained  if  a  significant  fraction  of  the  observed  subgrid  transfer  in  the 
physical  space  is  a  result  of  an  advection  of  the  small  scales  by  the  large  ones.  Even  though  such 
a  process  does  not  contribute  to  the  change  in  the  energy  of  the  large  scales  integrated  over  entire 
domain,  it  affects  their  energy  locally  in  the  physical  space  and  constitutes  a  legitimate  part  of  the 
transfer  in  the  physical  space  according  to  the  definition  (11).  Obviously,  such  a  process  will  be 
most  pronounced  at  the  locations  where  the  advective  velocity  is  large  and  these  locations  coincide 
with  the  regions  of  high  large  scale  energy. 

In  Fig.  6  the  subgrid-scale  transfer  is  overlayed  on  the  large  scale  dissipation  field  (15).  As  in 
the  case  of  correlations  with  the  enstrophy,  the  regions  of  the  large  transfer  are  generally  found 
on  peripheries  of  the  regions  with  large  dissipation.  This  is  perhaps  not  surprising  since  both  the 
enstrophy  and  the  dissipation  are  determined  by  the  strength  of  the  velocity  gradients.  However, 
it  should  also  be  noted  that  the  regions  of  the  large  enstrophy  and  the  large  dissipation  are  not  in 
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the  same  locations 


,  with  the  dissipation  largest  in  the  regions  of  small  enstrophy  vain 's.  and  vn  e 

versa. 

The  vorticity  and  the  enstrophy  production  terms,  (16)  and  (17),  respectively,  we  e  found  to  be 
fairly  well  correlated  with  the  enstrophy  field.  This  in  turn  implies  that  .lie  subgrid-scale  transfer 
will  be  largest  at  the  peripheries  of  the  regions  with  large  values  of  these  quantities,  similarly  to 
the  case  of  the  enstrophy.  Plots  of  the  transfer  superimposed  on  the  contour  plots  of  the  vorticity 
and  the  enstrophy  production  terms  are  shown  in  Figs.  7  and  8,  respectively,  and  confirm  this 
supposition.  Finally,  in  Fig.  9  the  comparison  is  made  between  the  subgi id-scale  transfer  and 
the  second  invariant  for  a  higher  resolution*5  held,  run  B.  Again,  we  find  that  the  correlations 
between  both  fields  is  such  that  the  most  intense  transfer  is  located  around  the  edges  of  regions 
with  significant  values  of  the  second  invariant. 

With  the  exception  of  the  energy,  all  the  above  quantities  are  various  combinations  of  the  first 
order  velocity  derivatives.  All  these  combinations  correlate  poorly  with  the  computed  subgrid- 
scale  energy  transfer  in  a  sense  that  large  values  of  these  quantities  are  found  in  different  spatial 
locations  than  large  positive/negative  values  of  the  transfer.  However,  the  transfer  is  usually  most 
intense  in  the  regions  located  on  peripheries  of  the  regions  characterized  by  large  values  of  these 
quantities.  It  may  be  expected  that  at  these  locations  gradients  of  the  above  analyzed  quantities 
will  be  large.  It  suggests  that  correlations  between  the  subgrid-scale  energy  transfer  and  quantities 
involving  second  order  derivatives  of  the  velocity  field  may  be  better  than  Lr  quantities  based  on 
the  first  derivatives.  One  possible  quantity  of  this  kind,  which  is  of  interest  in  the  context  of  this 
work,  is  the  classical  expression  of  Smagorinsky  for  the  subgrid-scale  ener^v  transfer.  In  Fig.  10 
we  plot  this  quantity  superimposed  on  the  exact  subgrid-scale  energy  tram  er.  The  Smagorinsky’s 


model  captures  locations  of  the  most  intense  transfer  much  better  than  the  other  invr.stigat<-«i 
quantities,  with  the  exception  of  the  energy.  The  obvious  drawback  of  this  model  is  that  it  has  a 
purely  dissipative  character  predicting  the  energy  transfer  always  from  the  large  to  the  small  scales 
whereas  the  exact  transfer  contains  regions  of  the  forward  as  well  as  inverse  transfer,  both  being 
significant. 

5  Conclusions 
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Figure  Captions 

Figure  1.  Isosurface  plot  of  the  total  enstrophy  for  the  velocity  field  A. 

Figure  2.  Cut  through  the  enstrophy  field  in  a  horizontal  plane:  (a)  located  at  z  =  Lj 4;  (b) 
located  at  z  =  L/ 2,  where  L  is  the  periodicity  length. 

Figure  3.  The  subgrid-scale  energy  transfer  in  a  horizontal  plane:  (a)  located  at  z  =  L/4;  (b) 
located  at  z  =  L/2,  where  L  is  the  periodicity  length. 

Figure  4.  The  subgrid-scale  energy  transfer  superimposed  on  the  large  scale  enstrophy  in  a 
horizontal  plane  z  =  Lj 4. 

Figure  5.  The  subgrid-scale  energy  transfer  superimposed  on  the  large  scale  energy  in  a  hori¬ 
zontal  plane:  (a)  z  =  L/4;  (b)  z  =  L/2. 

Figure  6.  The  subgrid-scale  energy  transfer  superimposed  on  the  large  scale  dissipation  in  a 
horizontal  plane  z  =  L/4. 

Figure  7.  The  subgrid-scale  energy  transfer  superimposed  on  the  large  scale  vorticity  production 
in  a  horizontal  plane  z  =  L/4. 

Figure  8.  The  subgrid-scale  energy  transfer  superimposed  on  the  large  scale  enstrophy  produc¬ 
tion  in  a  horizontal  plane  z  =  L/4. 

Figure  9.  The  subgrid-scale  energy  transfer  superimposed  on  the  second  invariant  of  the  large 
scale  velocity  gradient  tensor  (the  velocity  field  B). 

Figure  10.  The  exact  subgrid-scale  energy  transfer  superimposed  on  the  transfer  obtained  from 
the  Smagorinsky’s  eddy- viscosity  model  in  a  horizontal  plane  z  —  Lj 4. 
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